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Chapter 1 

Isolated impurity states 



1.1 Shallow Impurities 

A typical energy level diagram is shown on Fig. [TTT] Shallow levels allow a universal de- 
scription because the spread of wave function is large and the potential can be treated as 
from the point charge, 

U(r) = e 2 1 ' kt . 

To find impurity states one has to treat Schrodinger equation (SE) including periodic 
potential + Coulomb potential of the defect. 

Extremum at the center of BZ 

Then for small k we have 

We look for solution of the SE 

(Ho + U)if> = EiJj 



e d 



Figure 1.1: band diagram of a semiconductor 
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CHAPTER 1. ISOLATED IMPURITY STATES 



in the form 



ip = £ n /(k')0 n / k /(r) , 



n'k' 



where 0„'k'( r ) are Bloch states. By a usual procedure (multiplication by 0* k (r) and inte- 
gration over r) we get the equation 

[E n (k) - E}B n (k) + Y,U:i,B n ,(k) = 



n'k' 



jjnk 



Then, it is natural to assume that B(k) is nonvanishing only near the BZ center, and to 
replace central cell functions u by their values at k = 0. These function rapidly oscillate 
within the cell while the rest varies slowly. Then within each cell 

/ in Ko u n'o dr = 5 nn i 
J cell 

because Bloch functions are orthonormal. Thus, 

[E n (k) - E]B n (k) + J2U(kk')B n (k') = 

n' 

Aire 2 



U(kk) = y J e l ^- bk ' )r U(r) dr 



KVlk-k'l 2 



Finally we get 



2 U2 



2m 



-E 



k-k 



7^«( k ') 



where one can integrate over k in the infinite region (because B n (k) decays rapidly. 
Coming back to the real space and introducing 

1 



F(r) 



we come to the SE for a hydrogen atom, 



v Y.B n {k)e 



ikr 



a 2 n2 e 2 

_ V 2 

2m kt 



Here 



Et 
F(r) 



1 e 4 m 
T 2 2k 2 U 2 ' 



Fir) = EF(r) 



t = 1,2... 



[na 



3N-1/2 



exp(— r/a), a = h 2 n/me 2 



For the total wave function one can easily obtain 

tj} = u n0 {r)F(r) . 
The results are summarized in the table. 
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Material 


K 


m/rriQ 


Eu (th.) 
(meV) 


E Xs (exp.) 
(meV) 


GaAs 


12.5 


0.066 


5.67 


Ge:6.1 Si: 5.8 
Se: 5.9 S: 6.1 
S: 5.9 


InP 


12.6 


0.08 


6.8 


7.28 


CdTe 


10 


0.1 


13 


13.* 



Table 1.1: Characteristics of the impurity centers. 



Several equivalent extrema 

Let us consider silicon for example. The conduction band minimum is located at k Z Q = 
0.85(27r/a) in the [100] direction, the constant energy surfaces are ellipsoids of revolution 
around [100]. There must be 6 equivalent ellipsoids according to cubic symmetry. For a 
given ellipsoid, 




Here rri£ = 0.916 mo, m t = 0.19 mo- According to the effective mass theory, the energy lev- 
els are iV-fold degenerate, where n is the number of equivalent ellipsoids. In real situation, 
these levels are split due to short-range corrections to the potential. These corrections 
provide inter-extrema matrix elements. The results for an arbitrary ratio 7 = m t /m^ can 
be obtained only numerically by a variational method (Kohn and Luttinger). The trial 
function was chosen in the form 

F = (7ra||a^_) _1//2 exp 

and the parameters were chosen to minimize the energy at given 7. Excited states are 
calculated in a similar way. The energies are listed in table |1.1| . 



Material 


Ei, (meV) 


E 2po (meV) 


Si (theor.) 




31.27 




11.51 


Si(P) 


45.5 


33.9 


32.6 


11.45 


Si(As) 


53.7 


32.6 


31.2 


11.49 


Si(Sb) 


42.7 


32.9 


30.6 


11.52 


Ge(theor/) 




9.81 




4.74 


Ge(P) 


12.9 


9.9 




4.75 


Ge(As) 


14.17 


10.0 




4.75 


Ge(Sb) 


10.32 


10.0 




4.7 



Table 1.2: Donor ionization energies in Ge and Si. Experimental values are different 
because of chemical shift 



x 2 + y 2 z 2 

— H 

2 ' 
a , an 



1/2 
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CHAPTER 1. ISOLATED IMPURITY STATES 



Impurity levels near the point of degeneracy 

Degeneracy means that there are t > 1 functions, 

0i k , j = i,2..t 

which satisfy Schrodinger equation without an impurity. In this case (remember, k « 0), 

v = E i? i( r )^o( 6r )- 



The functions satisfy matrix equation, 



t 

E 

i'=i 



E H$PrtP + U(r)Sj 

a,P=l 



Fj< = EFj . 



If we want to include spin-orbital interaction we have to add 

1 



7~£so 



AtjiqC 2 



[sxVV]'P 



Here s is the spin operator while V is periodic potential. In general 7^-matrix is com- 
plicated. Here we use the opportunity to introduce a simplified (the so-called invariant) 
method based just upon the symmetry. 

For simplicity, let us start with the situation when spin-orbit interaction is very large, 
and split-off mode is very far. Then we have 4-fold degenerate system. Mathematically, it 
can be represented by a pseudo-spin 3/2 characterized by a pseudo- vector J. 

There are only 2 invariants quadratic in p, namely p 2 and (p • J) 2 . Thus we have only 
two independent parameters, and traditionally the Hamiltonian is written as 



n 



i 

m 



2 



7i 



+ ~^(P ' J ) 



;i.2) 



That would be OK for spherical symmetry, while for cubic symmetry one has one more 
invariant, Y^ipfJf- As a result, the Hamiltonian is traditionally expressed as 



n 



i 

m 



2 



7i + 



73(p-J) 2 + (73 - 7 2 )EA 2 ^ 



;i.s) 



This is the famous Luttinger Hamiltonian. Note that if the lattice has no inversion center 
there also linear in p terms. 

Now we left with 4 coupled Schrodinger equations (|1.1| ). To check the situation, let us 
first put U(r) — and look for solution in the form 



Fj = A,(k/A;)e lkr , k = |k 
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The corresponding matrix elements can be obtained by substitution Kk instead of the 
operator p into Luttinger Hamiltonian. The Hamiltonian (|1.2|) does not depend on the 
direction of k. Thus let us direct k along z axis and use representation with diagonal Jf . 
Thus the system is decoupled into 4 independent equation with two different eigenvalues, 



Ep 



2mn 



Et 



2m n 



If 



7i ± 2 7 > 

both energies are positive (here the energy is counted inside the valence band) and called 
the light and heavy holes. The effective masses are 

mp(h) = m / ( 7l ± 7) . 



The calculations for the full Luttinger Hamiltonian ( |1.3|) require explicit form of J-matrices. 
It solutions lead to the anisotropic dispersion law 



2m 



{^k 2 ± 4 [ 72 2 fc 4 



+12(7 3 2 - 7 2 2 )(^ + A;^ + A;^ 



1/2 



The parameters of Ge and Si are given in the Table [TTT] 



Material 


7i 


72 


73 


A 


K 


Ge 


4.22 


0.39 


1.44 


0.044 


11.4 


Si 


13.35 


4.25 


5.69 


0.29 


15.4 



Table 1.3: Parameters of the Luttinger Hamiltonian for Ge and Si 

The usual way to calculate acceptor states is variational calculation under the spherical 
model flOp. In this case the set of 4 differential equations can be reduced to a system of 
2 differential equation containing only 1 parameter, (3 = mi/rrih. 



Asymptotic behavior of the impurity- state wave function 

At large distances the wave functions of the localized states decrease exponentially, and for 
many problems it is sufficient to know only the exponential factor. Let us restrict ourselves 
by the effective massmethod which is valid at the distances much larger than the lattice 
spacing. For the case of non-degenerate band the Schrodinder equation for the envelope 
function, F, is of the form 

[H(-ihV) + U(r)}F = EF . (1.4) 
The asymptotic behavior can be understood within quasiclassical approximation, 



F = exp[iS(r)/K\ 
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where deep under the barriers the imaginary part of the action S is large.Q Using the 
above form for F and following the rules of quasiclassical approach where only the first 
derivatives of S should be kept we get a classical Hamilton- Yacoby equation 

H(VS) + U(r) =0. 

Now let us assume that we know both the action S(r) and its gradient VS* at a surface 
a. Now we can draw a ray through each point of a which satisfies classical equations of 




Figure 1.2: A quasi-classical trajectory passing through the point . 
motion and has p = VS* on the surface a. Then at arbitrary point r we have 

5(r) = 5(r CT )+ f p • dr' . 
According to classical mechanics, we have the Hamilton equations 

*= — P = -— • (1.5) 
dp ' dr 

Since by definition S^r^) <C S(r), and outside the range of the potential p = p = const, 

S(r) = p • r 

The group velocity v = dE/dp corresponding to the vector p should be perpendicular to 
the surface a at the incident point, while its absolute value is determined by the condition 

H(p ) = E. (1.6) 

As a result, at large distances the wave function is a plane wave 

p _ e «po-r/?i _ 

1 The condition for the validity of the quasiclassical approach reads 
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Consider the simplest case, 7Y(p) = p 2 /2m where the group velocity v = V p 7i || p. Then 
for negative energies we can write 



>2m\E\ 

F(r) = e- qr , q 



h 

This solution appears exact for the hydrogen-like spectrum with 

q = a -1 , a = h 2 K/me 2 . 

The quantity a is the effective Bohr radius. Remember that for the validity of the approx- 
imation under consideration the inequality 

qa < 1, a = h 2 /m e 2 

must be met. 

Now let us consider the case of ellipsoidal spectrum, 

W(p) = *Lt*l + A . (1.7) 
2mi irrit 



Introducing the mass ratio 7 = m t /mi and the quantity d = ^p\ x + Po y + l 2 p\ z we have 
the set of equations 

Pox = n x d, poy = n y d,p 0z = n z d/j , (1.8) 
Substituting these relations to Eq. (|1.6|) for negative E we get 

d 2 _ 2m t \E\ 

n 2 x + n 2 y + n 2 J~f ' 

As a result, 



g(n) = /r\/2mi|£|(n2 +n 2 v + n 2 z h) . (1.9) 



We observe that the shape of the wave function is also ellipsoidal. In particular, for Ge we 
have qf 1 = 13.8 A (n x = n y = 0), q^ 1 = 61.4 k{n z = 0). 

As we have mentioned, in semiconductors with several ellipsoids (Ge, Si) the wave 
function is represented by a linear combination of the functions corresponding to different 
ellipsoids with different ^(n). At very large distances we are interested in 

q(n) = min{(7j(n)}. 

We do not discuss here the case of degenerate bands. One can find details in the book 
by Shklovskii and Efros [□] 
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Chapter 2 



Localization of electronic states 



What happens when the number of doping impurities is large and their wave functions 
overlap? 

2.1 Narrow bands and Mott transition 

As a simple example, consider an ordered lattice of impurities, the potential being 

V(r) = J2U(r-r 3 ). 

j 

Assume that we know the eigenstates <p n and eigenvalues, E n , of a single-impurity problem. 
Here we shall use single-band approximation, and therefore ignore Bloch (central-cell) 
factors. Also, we assume the impurity bandwidth is less that the spacing between E n and 
restrict ourselves by the lowest level, E . Along the tight-binding method, we can expand 
the wave functions in terms of the above eigenvalues, 

*jj = J2a j (j)(r-r j ), £ |a/ = 1 , 

j j 

the energy being 

E = J2e ikm I(m) . 

m 

The only important property of the overlap integrals, I(m) is that it is small (tight-binding 
approximation! ) . 

In this way we get energy bands. For example, for a simple cubic lattice 

£ = 2/(& )£cos(fc i & ) 

i 

where here b denotes the lattice constant for impurity lattice. Expanding near k = we 
get 

E(k) = 6/(6 ) - I(bo)k 2 b 2 
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CHAPTER 2. LOCALIZATION OF ELECTRONIC STATES 



the effective mass being m* = h 2 /2Ib\. As the lattice constant increases, I (bo) oc exp(— flbo/a) 
where (3 is a number of the order 1. 

According to our model, each impurity adds one electron, and each state possesses 
2-fold spin degeneracy. Thus the impurity lattice is a metal. 

Is this correct? In principle, no, because we neglected electron-electron interaction. 
In fact, two electrons at the same site repel each other the energy can be estimated as 
Uq ~ e 2 /a. It Uq is comparable with the bandwidth ~ I (bo), then one must allow for the 
interaction. At large bo the band is narrow and this is the case. Let us plot the electron 
terms a a function of 1/bo The insulator-to-metal transition of this kind is usually called 



* E 




A l/b 



Figure 2.1: Dependence of electron bands on impurity sublattice period bo- To the left of 
point A is an insulator, to the right a metal. 

the Mott transition. 



2.2 Anderson transition 

We come back to single-electron approximation and concentrate on disorder. Suppose that 
impurities are still ordered, but depths of the potential wells are random, Fig. (|2.2|) The 
effective Hamiltonian has the form 

H = Y^ e j a t a J + H J ( m ) a f a j+m ■ 

The energies e are assumed to be uncorrelated and random, the distribution being 

P(A _ [ 1/A , \e\ < A/2 
1 1 ~ 1 , |e| > A/2 

Anderson has formulated the following criterion for localization. Let us place a particle at 
the site i and study its decay in time due to quantum smearing of the wave packet. The 
state is called extended if \^i(t)\ 2 — > at t — » oo. 
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Figure 2.2: Potential wells in the Anderson model. 

The Anderson model has one dimensionless parameter, A/ 1, where / is the overlap 
integral for nearest neighbors. In a 3D system there is a critical value, A c /I, at which 
delocalized states begin to appear in the middle band. In ID case the states are localized 
at any disorder. 2D case is a marginal, the states being also localized at any disorder. 

2.2.1 Two-state model 

Now let us turn to 3D case and try to discuss Anderson result in a simplified form. Namely, 
consider two potential wells. When isolated, they can be characterized by the energies e\ 
and e 2 and wave functions 4>i an d <t>2- At e\ = e 2 the wave functions are 



= ^(01 ±02), 

the energy difference being Ej — Eji = 21. This solution is reasonably good at |ei — e 2 | <C /. 
In general, 

$1,11 = Cl01 ± c 2 2 , 

and we have a matrix Schrodinger equation, 

( A/2-E I \ 

\ I* -A/2-E ) ' 

Here the origin for energy is taken at (ei + e 2 )/2, while A = e\ — e 2 . The secular equation 
is thus 

E 2 - (A/2) 2 - \I\ 2 = -> Ejjj = ± v /(A/2) 2 + |/| 2 . 

Consequently, 

Ei — En = V / A 2 + 4|/| 2 . (2.1) 

The ratio 



C2 A/2± J(A/2) 2 + / 2 
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Thus at A ^> / either c\ or C2 is close to 1, and collectivization does not occur. 
Now, following Thouless, let us chose a band 

5/2 <e< 5/2, 5~I 

and call the sites resonant if their energy falls into the band. Then let us look for connected 
resonant states which share a site. Non-resonant sites can be disregarded as I <C A. 

It is clear that it must be a threshold in the quantity A/ 1 where the transition takes 
place. If one assumes that the connected cluster is a ID chain, then the bandwidth is 47. 
In such a model, 

Ac = ± 

I X c 

where x c is the percolation threshold for the site problem. This is quite a good estimate, 
see the Table |2J] Finally, one arrives at the following profile of density-of-states (DOS) 



Lattice 


x c 


4/x c 


A c /I 


Diamond 
Simple cubic 


0.43 
0.31 


9.3 
12.9 


8.0 
15 



Table 2.1: Percolation thresholds and critical values A c j I obtained numerically for different 
lattices. 

for the Anderson model, Fig. According to this picture, there both extended and 





g(E) 


-E c 


E 
c 


E 



Figure 2.3: Density of states in the Anderson model. Localized states are shaded. Mobility 
edges are denoted as E c . 

localized states separated by the mobility edges. If we vary the electron density the Fermi 
level will move with respect to the bands and it may cross the mobility edge. At this 
point a transition in conductance takes place which is conventionally called the Anderson 
transition. 

To illustrate the Anderson transition let us discuss the case of an amorphous semi- 
conductor. The presence of a short-range order preserves the band picture, so there are 
forbidden gaps which separate the allowed regions of energy. However, structural "faults" 
give rise to the tail in the DOS, and the band boundaries appear smeared. If the Fermi 
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Ey H E c E 



Figure 2.4: DOS in an amorphous semiconductor. Localized states are shadowed. 

level falls into the localized region the transport is possible only due to thermal activation 
into the delocalized states region, or by hopping between the localized states. Both mecha- 
nisms lead to an exponential temperature dependence of conductance at low temperatures. 
At very low temperatures a weak tunneling between the localized states is the only source 
of transport, so it becomes essentially temperature-independent. 

If we increase the electron concentration, the Fermi level moves to the direction of delo- 
calized states, and the insulator-to metal-transition takes place. As a result, a temperature- 
activated hopping crosses over to a weak temperature dependence which is called the metal- 
lic conductance. Such transition is called the Mott-Anderson transition because it bears 
the features of both models. 

Anderson localization appears important for metal-insulator-semiconductor (MIS) struc- 
tures which are extensively used in modern electronics. Such a structure is is most often 
realized with a SiC>2 film (insulator) sandwiched between Si substrate and a flat metal 
electrode. Electrostatic forces force the energy bands to bend to keep a constant electro- 
chemical potential. As a result, the changes re-distribute creating on the semiconductor 
surface a narrow inversion layer, i. e. the layer with the carriers of the type opposite to 
the ones in the bulk, see Fig. |2.5| . By changing the applied potential one can control con- 



Figure 2.5: Band diagram near the MIS interface. The inversion layer is shadowed, fj, is 
the chemical potential. 

ductivity of the inversion layer over a wide range. This principle is used for the so-called 
field-effect transistors (FET). At low carrier concentrations the conductance of the layer is 
activated. As the concentration increases, the Fermi level fj, moves towards the extended 
states, the conductivity increases and becomes a weak function of the temperature. 




b 



d 
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Now a question arises: Is the Mott- Anderson transition abrupt or continuous? The 
first answer to this question has been suggested my Mott who arrived at the concept of 
minimal metallic conductivity. This concept is close to the so-called Ioffe-Regel principle 
that the mean free path I which enters the Drude formula for 3 dimensional electron gas 
(3DEG), 

_ e 2 p 2 F £ 

cannot be less than the de Broglie electron wavelength, H/pp. In this way we obtain 

e 2 

u mm f-n 

Within the Anderson model, £ mm = a$ where ao is the lattice constant, while for the Mott 
model £ min corresponds to the typical distance between the impurities. According to Mott's 
suggestion, at T — > the conductivity reaches its minimal value, and then drops to zero, 
see Fig. |2.6| . Mott's arguments for the two-dimensional case lead to the estimate 




Figure 2.6: Mott's suggestion of the metal-to-insulator transition. Zero-temperature con- 
ductivity as a function of the Fermi level \i. 

e 2 
n 

The modern theory contradicts to the Mott's concept, it predicts a continuous transition. 



2.3 Modern theory of Anderson localization 

We start from the case of metallic conductance to find in which way disorder influences 
the conductance. Then we shall discuss the aspects of strong localization. 



2.3.1 Weak localization 

Consider noninteracting electrons having ppl 3> h and passing between the points A and 
B through scattering media. The probability is 

2 



w 



i i^j 



(2.2) 
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Here Ai is the propagation amplitude along the path i. The first item - classical probability, 
the second one - interference term. 

For the majority of the trajectories the phase gain, 

rB 

A(p = h~ 1 pd\ > 1 , (2.3) 

J A 

and interference term vanishes. Special case - trajectories with self-crossings. For these 




Figure 2.7: Feynman paths responsible for weak localization 

parts, if we change the direction of propagation, p — > — p , dly — > —dl, the phase gains are 
the same, and 

\Ai + A 2 \ 2 = |A| 2 + \A 2 \ 2 + 2A ± A* = 4|Ax| 2 . 

Thus quantum effects double the result. As a result, the total scattering probability at 
the scatterer at the site O increases. As a result, conductance decreases - predecessor of 
localization. 

Probability of self-crossing. The cross-section of the site O is in fact de Broglie electron 

wave length, A ~ h/p F . Moving diffusively, if covers the distance \fxf ~ VTJi, covering the 
volume (Dt) d / 2 b 3 ~ d . Here d is the dimensionality of the system, while b is the "thickness" 
of the sample. To experience interference, the electron must enter the interference volume, 
the probability being 

v\ 2 dt 
( Dt y/2 b 3-d ■ 

Thus the relative correction is 




Figure 2.8: On the calculation of the probability of self-crossing. 



18 



CHAPTER 2. LOCALIZATION OF ELECTRONIC STATES 



The upper limit is the so-called phase-breaking time, t v . Physical meaning - it is the 
time during which the electron remains coherent. For example, any inelastic or spin-flip 
scattering gives rise to phase breaking. 

3d case 

In a 3d case, 



Ad v\ 2 ( 1 1 \ /A\ 2 A 2 



a D 3 / 2 \y/r J^J \£J £L, 



(2.5) 



Here we have used the notations 



D~vt, t~£/v, L v = ^Dr v . (2.6) 
On the other hand, one can rewrite the Drude formula as 

n e e 2 r n e e 2 £ e 2 p 2 F £ 



a 



m p F h 

In this way, one obtains 

Important point: The second item, though small, is of the most interest. Indeed, it 
is temperature-dependent because of inelastic scattering. There are several important 
contributions: 

• Electron-electron scattering: 

~ he F /T 2 . 

Thus 

a(T)-a(0) f\\ 3/2 {T 



<r(0) \£ J \e F , 

It is interesting that at low temperatures this quantum correction can exceed the 
classical temperature-dependent contribution of e-e-scattering which is oc T 2 . It is 
also important to note that the above estimate is obtained for a clean system. It 
should be revised for disordered systems where electron-electron interaction appears 
more important (see later). 



Electron-phonon-interaction. In this case,Q 



D 



and 

a(T)-a(0) _ f\\ 3/2 f T \ 1/2 ( T 



cr(0) \£ J \e F J \%L0d 

1 Under some limitations. 
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Thus there is a cross-over in the temperature dependencies. To obtain the dominating 
contribution one has to compare r" 1 . Consequently, at low temperatures e-e-interaction 
is more important than the e-ph one, the crossover temperature being 

T ~ (huip/ep) ~ IK. 



Low-dimensional case 

If the thickness b of the sample is small 0, then the interference volume X 2 v dt has to be 
related to (Dt) d ^ 2 b 3 ~ d . For a film d = 2, while for a quantum wire d = 1. As a result, 

Aa ^ _v\ 2 j fe-Mn^/r), d = 2 



D 1 d = 1 



It is convenient to introduce conductance as 

G = ab 3 - d . 

We have, 

h [ L v , d = 1 

Important note: In low- dimensional systems the main mechanism of the phase breaking 
is different. It is so-called low-momentum-transfer electron-electron interaction which we 
do not discuss in detail 



2.3.2 Weak localization in external magnetic field 

In a magnetic field one has to replace p — > p + (e/c)A (remember - we denote electronic 
charge as — e). Thus the product AA* acquires an additional phase 

2e f 2e f $ 

Am H = — A AdL= — i curl A dS = 4vr— (2.10) 

dh J chJ $ 

where $ is the magnetic flux through the trajectory while $o — 2irtic/e is the so-called 
magnetic flux quantum. Note that it is 2 times greater than the quantity used in the theory 
of superconductivity for the flux quantum. 

Thus magnetic field behaves as an additional dephasing mechanism, it "switches off" 
the localization correction, increases the conductance. In other words, we observe negative 
magnetoresistance which is very unusual. 

To make estimates, introduce the typical dephasing time, to get Aipn ~ 27T for the 
trajectory with L ~ \fDtn. In this way, 

t H ~%/(H D)~l 2 H /D, l H = ^ch/eH. (2.11) 

2 The general from of the criterion depends on the relationship between the film thickness, b, and the 
elastic mean free path £. The result presented is correct at b <C I. At L 3> b 3> t one can replace the lower 
limit t of the integral ( fO| ) by t& = b 2 /D. 
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The role of magnetic field is important at 

t H < TV -> H>H ~ $o/(Dt v ) « ^c/Lj . 

Note that at H ^ H 

uj c t ~ h/(e F r v ) < 1 

that means the absence of classical magnetoresistance. Quantum effects manifest them- 
selves in extremely weak magnetic fields. 

To get quantitative estimates one has to think more carefully about the geometry of 
diffusive walks. Let consider the channel of 2DEG with width W. The estimates given 
above are valid for £, L v <C W, exact formulas could be found, e. g. in Ref. [□]. Usually H 
is very weak, at L v = 1 /zm Hq m 1 mT. The suppression of weak localization is complete 
at H > h/e£ 2 , still under conditions of classically weak magnetic field, uj c t ^ 1. 

The situation is a bit different at W < L v , this case can be mentioned as one- 
dimensional for the problem in question, see Fig. |2.9| . Then a typical enclosed area is 




Figure 2.9: On quasi ID weak localization. 

S ~ WLp, and the unit phase shift takes place at 

t H ~ L%/DW 2 , -> i/ ~ fe/e^L v . 
Some experimental results on magnetoresistance of wide and narrow channels are shown 

in Fig. Eng . 

There are also several specific effects in the interference corrections: 

• anisotropy of the effect with respect of the direction of magnetic field (in low- 
dimensional cases); 

• spin-flip scattering acts ClS db dephasing time; 

• oscillations of the longitudinal conductance of a hollow cylinder as a function of 
magnetic flux. The reason - typical size of the trajectories. 
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8(G) 



Figure 2.10: Experimental results on magnetoresistance due to 2D weak localization (upper 
panel) and due to ID weak localization in a narrow channel (lower panel) at different 
temperatures. Solid curves are fits based on theoretical results. From K. K. Choi et ai, 
Phys. Rev. B. 36, 7751 (1987). 



2.3.3 Aharonov-Bohm effect 

Magnetoconductance corrections are usually aperiodic in magnetic field because the in- 
terfering paths includes a continuous range of magnetic flux values. A ring geometry, 
in contrast, encloses a continuous well-defined flux $ and thus imposes a fundamental 
periodicity, 

= G($ + n$ ) , $ = 2nhc/e . 

Such a periodicity is a consequence of gauge invariance, as in the originally thought exper- 
iment by Aharonov and Bohm (1959). The fundamental periodicity 

eb 

comes from interference of the trajectories which make one half revolution along the ring, 
see Fig. [2.1 1| . There is an important difference between hc/e and hc/2e oscillations. The 



first ones are sample-dependent and have random phases. So if the sample has many 
rings in series or in parallel, then the effect is mostly averaged out. Contrary, the second 
oscillations originate from time-reversed trajectories. The proper contribution leads to a 
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Figure 2.11: Illustration of Aharonov-Bohm effect in a ring geometry, (a) Trajectories 
responsible for hc/e periodicity, (b) trajectories of the pair of time-reversed states leading 
to /ic/2e-periodicity. 



minimum conductance at H = 0, thus the oscillations have the same phase. This is why 
/ic/2e-oscillations survive in long hollow cylinders. Their origin is a periodic modulation of 
the weak localization effect due to coherent backscattering. Aharonov-Bohm oscillations 
in long hollow cylinders, Fig. |2.12| , were predicted by Altshuler, Aronov and Spivak [Ef 
and observed experimentally by Sharvin and Sharvin [0]]. A rather simple estimate for the 
magnitude of those oscillations can be found in the paper by Khmelnitskii [0]. Up to now 




Figure 2.12: On the Aharonov-Bohm oscillations in a long hollow cylinder. 



we have analyzed the localization corrections ignoring electron-electron interaction. In the 
following section we shall discuss the influence of disorder on e — e interaction. 



2.3.4 Electron-electron interaction in a weakly disordered regime 

Let us discuss the effect of the e — e interaction on the density of states. Let us concentrate 
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Direct Exchange 



Figure 2.13: On the calculation of e — e interaction. 



on the exchange interaction, shown in a right panel of Fig. [2.13| 

d 3 k 



Ae 



|p-^k|< PF (2-nhy 



(2.12) 



Here g(k) is the Fourier component of the interaction potential, the sign "-" is due to 
exchange character of the interaction. In the absence of screening g(k) = 4ne 2 /k 2 , and 



Ae 



-Aire 



-A-ne 2 Ti 2 

„2 / 



, o d k 

P -hk\< PF (2tt) 3 

1 d 3 p' 



Pf 



<pf |p - P'| 2 (27rfi) J 
P 2 -P 2 f, P + Pf^ 



2p 



In 



P~PF, 



To obtain this formula it is convenient to use spherical coordinates - (dp') 
and auxiliary integral 



d(cos0) 



i p 2 + p' 2 — 2pp' cos i 



ppi 



In 



p + p' 
p — p' 



(2.13) 
2irp' 2 dp' d(cos 4>) 

(2.14) 



At small p — pf it is convenient to introduce £ = v (p — pp) "C €f to get (omitting the 
irrelevant constant) 

Ae = -(e 2 /nhv)£\n(2p F v/0 . (2.15) 
Screening can be allowed for by the replacement k~ 2 — > (k 2 + k 2 ) -1 that replaces p±p' in 



(2.16) 



the argument of the logarithm in Eq. (2.14) by y (p ±p') 2 + (Tin) 2 . As a result, 

Ae = — (e 2 /%hv )£ ln(2pF/hK) at TiK<^pp- 



This is a simple-minded estimate because it ignores the interference of the states. Indeed, 
if the two states differ in the energy by |£| the coherence time is Ti/\^\. If the electron 
returns back for this time, then the effective interaction constant increases by 
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Thus 



<ff = e 2 (l + a € ). 



In a similar way 



Av Ag 



v 



9 




Here g is the density of electronic states. Consequently, since a oc v, we arrive to a specific 
correction to conductance. Comparing this correction with the interference one we conclude 
that interaction dominates in 3d case, has the same order in 2d case and not important 
in Id case. Another important feature is that the interaction effects are limited by the 
coherence time h/\£\ ~ tt = h/ksT rather than by t v . Usually, t v ^> tt- Consequently, 
the interference effects can be destroyed by weaker magnetic fields than the interaction 
ones (important for separation of the effects). 

2.3.5 Scaling theory of Anderson localization 

If we come back to the interference correction for d — 2, 1 we observe that it increases with 
Tip, or at T — > 0. Thus the corrections becomes not small. We can also prepare the samples 
with different values of the mean free path t. 

What happens if the corrections are not small? Anderson (1958) suggested localization 
of electronic states at T = 0. This suggestion has been later proved for an infinite Id 
system, as well as for an infinite wire of finite thickness (Thouless, 1977). Later it has been 
shown that 



where L loc ~ £ in the first case and {bppjKf'i for the second one (exponential localization). 
It seems that such a law is also the case for a metallic film (rigorous proof is absent). 
Very simple-minded explanation - over-barrier reflection + interference of incoming 
and reflected waves. Because of the interference the condition T = is crucial (no phase 
breaking). This explanation is good for one-dimensional case. 

A little bit more scientific discussion. Consider interference corrections to the con- 
ductance at T — » 0. In this case one has to replace 



One can conclude that in 3d case the relative correction is ~ (h/ppi) 2 <C 1 (usually). 
However, at d — 1, 2 it increases with the size. 
At what size Aa m a? 



G oc exp(-L/L loc ) 



L 2 /D, L V ^L. 




(2.18) 



We can conclude that in 3d case localization takes place at ppi ~ h, while in Id and 2d 
case it takes place at any concentration of impurities. 
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Scaling hypothesis - According to the "gang of 4" (Abrahams, Anderson, Licciardello 
and Ramakrishnan) 

G(qL) = f[q,G(L)]. (2.19) 
Assuming q — 1 + a, a C 1, we can iterate this equation in a: 

G{L) = /[1,G(L)] ; 
aLG'(L) = a{df/dq) q=1 . (2.20) 

Denoting 

G-\df/dq) q=1 =P(G) 
we re- write the scaling assumption through the Gell-Mann & Low function, /3(G): 

dhiG/d\nL = (3{G). 

At very large G we can expect that the usual theory is valid: 

f S\/L = L, d = 3; 
G = a\ bL ± /L = b, d = 2; (2.21) 
[ b 2 /L, d = 1 

Thus, in the zero approximation, 

(3(G) = d- 2. 

Then we can use the weak localization approximation to find next corrections. One can 
get 

P(G)nd-2-a d G /G, (2.22) 

where 

G = e 2 /(ix 2 h) , a d ~ 1. 

Indeed, for d = 3 

lnG = ln[(a + Aa)L] ^\naL + (Aa)/a w In [(a + Aa\ L=00 ) L] + h 2 /(p 2 F £L) . 

Thus 



<91nG h ha G 

m= dlnL =1 -p^L = 1 'pjUG = 1 - a3 cY 
At small G one can suggest exponential localization: 



G ~ G exp(-L/L c ) -> /3(G) ~ ln(G/G ) . 

Thus we have the scenario shown in Fig. [2.14| . Believing in such a scenario we get local- 
ization for Id case (/3(G) < - conductance increases with the length). At d = 3 we have 
a fixed point at G c , which is unstable ((3 changes sign). Under the simplest assumptions 

/3(G) « 7 (lnG-lnG c )« 7 (G/G c -l), 
G = G*- - 1 at L = Lq (initial condition) 
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Unstable fixed point 

\. 1 




P(G) 


d=3 

G increases with L 


G decreases with L V, 


G d=2 1,1(1 


d=l 



Figure 2.14: Flow curves for Anderson localization. 
(Gr(°) is close to G c ) we obtain 



G 



(2.23) 



From Eq. (|2.23| ) one can find important dependencies near the critical point. It is natural 
to chose L « £ and to suggest that at this size 



cr ~ e p F £/h 



G^ = a L = (e 2 /h)(p F £/h) 2 . 



Now we can assume that we can control some parameter (say, impurity content, x, which 
effects the mean free path £), and that G^ is regular in this parameter. Then, at small x 
= G c (l + x). At L = £ we obtain 

G = G C {1 + x) {L/i)1 w G c exp[x(L/£) 7 ] . 

Of course, such an assumption valid only at 

x < 1 , L/£ > 1 . 

Thus at x < we obtain exponential localization with the characteristic length 



" 1/7 , x < . 



However, at x > G grows with L. In a spirit of the concept of phase transitions we can 
treat the quantity L c = £x~ 1 ^ 1 as a correlation length. At the scales of the order of L c the 
properties of conducting and insulator phases are similar. The above law can be valid only 
in the vicinity of G^ ~ G c . Then we match the usual Ohm's law, G = o~L. Consequently, 
the conductivity can be estimated as 

G(0) 1/7 

a ~ — - — ~ (TqX 1 , x > (J . 
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Thus, we predict power law. Experiments on the dielectric constant (k oc L 2 c ) support the 
value 



Don't forget that we discuss the case T = 0, and the conductance is supposed at zero 
frequency. The range of applicability is given by the inequalities 



This is almost impossible to meet these conditions, so usually people extrapolate experi- 
mental curves to T = 0, u = 0. It is a very subtle point because, as it was shown, the 
conductance at L < L^ is not a self- averaging quantity with respect to an ensemble of 
samples. More precise, the fluctuations between the samples cannot be described by the 
Gaussian law, their distribution being much wider. Then, 

• have all these scaling assumptions any sense? 

• Why they reasonably agree with the experiment? 

The answer is positive because both the scaling predictions and the experiment are valid 
as an extrapolation from the region L > L v . 

As the temperature grows, fluctuations decrease and the conductance tends to the 
Ohm's law. 

Role of e — e interaction 

Nobody can consider both disorder and interaction acting together. To get some under- 
standing of the role of e — e interaction let us consider a clean metallic conductor. Assume 
that under some external perturbation (say, pressure) the band overlap changes. In this 
way we control the Fermi level (number of electrons and holes). One can consider them 
as free ones as their kinetic energy p 2 F /2m* exceeds the potential energy e 2 /(/t f). In this 
way we come to the condition 



Otherwise electrons and hole form complexes - Wannier-Mott excitons. This state is insu- 
lating because excitons are neutral. This is only one of possible scenario. In general, the 
problem is a front end of modern condensed matter physics. 



7 = 0.6 ±0.1. 
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Chapter 3 

Impurity band for lightly doped 
semiconductors. 



The material is called lightly doped if there is only a small overlap between electronic states 
belonging to different impurities, Na 3 <C f . Here N is the impurity concentration while a 
is the localization length. In lightly doped materials conductivity vanished exponentially 
as T -> 0. 

Let us start with a material with only one type of impurities, say n-type. At T = each 
of donors must have an electron, the missing electron represents an elementary excitation. 
An excitation can be localized at any donor, the energies being slightly different. The 
reasons are 

• If one removes an electron, the remaining positive charges polarize the neutral donor 
located in vicinity. That contributes to donor ionization energy, the contribution 
being dependent on the configuration of neutral neighbors. 

• There is a quantum overlap between the donors being excited. 

The first mechanism is usually more important for realistic semiconductors having com- 
pensating impurities. 

Now let us assume that there are also acceptors which at capture electrons from the 
donors and are fully occupied T = 0. Thus we have neutral donors, negatively charged 
acceptors and an equal number of positively charged donors. Theses randomly distributed 
charges create a fluctuating random potential which localizes the electronic states. 

Let us count the energy e; of i-th donor from the energy of an isolated donor, —E . We 
have 

o2 



e 

K 



acc i don i _ 

I I 1 "* ~~ r 'l k^i \ Ti ~ rfc l 



The occupation numbers have to be determined to minimize the free energy (at T = - 
electrostatic energy). 

A typical dependence of the Fermi level /i on the degree of compensation, K = Na/N d , 
is shown in Fig. ^| Below we discuss limiting cases of small and large K. 
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Figure 3.1: Position of the Fermi level /i as a function of the degree of compensation, K. 

3.1 Low degree of compensation 

At K <C 1 most of donors keep their electrons, and only a small number is empty and 
charged positive. Each acceptor carries a negative charge, the concentration of charged 
donors is equal to Na- 

A positive charged hole will be located as close as possible to a negative acceptor. 
Thus, it occupies the donor the most closed to an acceptor. In a first approximation, each 
acceptor can be regarded as immersed in an infinite sea of donors. Transporting a hole 
from a donor situated near a charged acceptor to infinity requires work of order e 2 /nro, 
where is the average distance between the donors. 

To get a quantitative theory one has to add more input. Namely, there are some 
acceptor without holes, and some acceptors having 2 holes. Indeed, some acceptors do not 
have a close donors (Poisson distribution), so they bind holes very weakly. So holes prefer 
to find another acceptor with more close donors even at the expense to become a second 
hole. To illustrate the situation let us consider a configuration shown in Fig. |3.1| ,a. The 
work necessary to remove a hole equals to 

e 2 / kt — e 2 /2nr = e 2 /2nr , 



it becomes large at small r. 

It is curious that an acceptor cannot bind more than 2 holes. Consider a configura- 
tion shown in Fig. 3.1,b. The energy of attraction to the acceptor equals \f?>e 2 / ' nl while 
the repulsion energy from two other holes is 2e 2 /nl. We end at the following situation. 
There are 3 configurations - 0-complexes (negative) where there is no ionized donor near a 
particular acceptor, 1-complexes (neutral), and 2-complexes (positive). So the neutrality 
condition which fixes the Fermi level is actually 



No(fM) = NM 
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Figure 3.2: Donor-acceptor configurations 



For 0-complex, there must be no donors within the sphere of radius r M = e 2 / k/i from a 
fixed acceptor, the probability being exp(— AnN D r 3 /?>). Thus, 

It is much more difficult to find number of 2-complexes. Let us estimate it from above, 
N^in), as a concentration of pairs of donors whose additional energies e\ and e 2 exceed /x 
when 6ot/i donors are ionized. This estimate is larger because 

• it could be another close donor which initiates a 1-complex, 

• it can be more than one pair of donors near a given acceptor. 

Let us put the origin of coordinates at acceptor site and suppose that a donor is located 
at ri. The average number of donors in the element dr 2 is equal to N D dr 2 . Thus, the 
number of pairs with r 2 > r\ is 



N D j dr 2 9[ei(ri, r 2 ) - /i]0[e 2 (r 1 , r 2 ) - fj] . 



Here 

e 2 

Cl,2 = — 

k L|ri )2 | |ri - r 2 | 

To get the total concentration of pairs one has to multiply it by N D and integrate over r 1; 
and finally multiply by Na- As a result, 

N>(n) = N A N 2 D [ tfri / dr 2 e[ei(r!, r 2 ) - //]0[e 2 ( ri , r 2 ) - fj] . 

This estimate is very close to exact result. Integration yields 

N>(fi) = 7.14 • lO~ 4 N A (A7cN D rl/3) 2 . 

Using this estimate and solving the neutrality equation one obtains 

li = 0.61e D , e D = e 2 /nr D = (e 2 //t)(47riV D /3) 1/3 . 

Now let us discuss the shape of the peak. The majority of donors are far from few 
acceptors, so there must be a sharp peak at —E . The tails of g(e) above the Fermi level 
falls with the characteristic energy eo- Consequently, near the Fermi level 

g(e) oc N A /e D . 
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Long-range potential 

Above we did not take into account electrostatic interaction between the complexes. It is 
small indeed because N Q <C Njj. However, the interaction leads to an additional dispersion 
of the peak. 

An average fluctuations of charge produce an average potential of order e 2 N^ 3 /n. 
We shall show that long-range fluctuations are much more important. Let us introduce 
fluctuations of the complex concentrations for 0- and 2-complexes, 

&(r) = JVi(r) - (JVi(r)) . 

We consider the fluctuations as uncorrelated, 

(£o(r)£o(r')> = (6(r)6(r')> = N 8(v - r'), 
<£o(r)6(r')> = 0. 

Consequently, the charge fluctuations are 

(p(r)p(r / )) = 2e 2 iV 5(r-r / ). (3.1) 

Let us now consider a sphere of radius R where there is ~ N Q R 3 complexes. The 
typical charge in the sphere is e 2 (N R 3 ) 1 ^ 2 , the resultant potential being e 2 (N R 3 ) 1 ^ 2 / R. 
It diverges as R — > oo, and one has to allow for screening. We will check later that the 
screening potential varies little over a distance between the complexes. Therefore, one can 
employ the self-consistent field approximation. 

As a result, we and at the Poisson equation 

Att Awe 
= p(r) [N 2 (fi + e4>) -N (jm + e<f>)] . 

K K 

Here p(r) is a Gaussian random function whose correlator is given by (|3.1|) . This approach 
is consistent at 

|e0(r)| <C n 

when only small number of complexes are responsible for screening. At the same time, the 
majority of complexes contribute to charge fluctuations and are uncorrelated. In such a 
way we obtain the equation 

A0 = -^ + 4 (3.2) 

where the screening length ro is given by the expression 

1 Ane d 

^ = — — [N 2 {fi + e<P) - N (fi + e0]^ =o . 

Substituting the expressions for the complex concentrations we obtain 

r = O.SSA^ 172 ^ 76 = 0.58iV D 1/3 fr 1/2 . 
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Then one can calculate the potential distribution from Eq. (|3.2|) as 



where 



We have 



e0 r 



J p(r') K{r - r) dr' 



K(v) 



-r/r 



(e0(r)) 2 = | dn J dr 2 (p(r 1 )p(r 2 ))K(r - r x )K(r - r 2 ) 



2e 2 iVn / K 2 (r) 



Kr 



27rAT or 3 . 



(3.3) 



(3.4) 



(3.5) 



In a standard way, one can determine the distribution function near the Fermi level. It is 
Gaussian, 



7V 71 " 



exp[-(e0) 2 / 7 2 ] , 7 = 2e 2 (0(r)> 2 = O^K 1 / 4 . 



All the above results are applicable only at very small K. Indeed, to get a Gaussian 
fluctuations one has to assume r ^> Nq 1 ^ 3 or 

477^0^/3 > 1 . 

The condition can be rewritten as K^ 1 ^ 2 <C 0.05. 

A typical sketch of the impurity band for a weakly compensated semiconductor is shown 
in Fig. |3T3| 

c-band 



-e- -e-e- 




g(E) 



v-band 

Figure 3.3: Energy diagram for a weakly compensated semiconductor 
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3.2 High degrees of compensation 



Now we turn to the case 



1-K<1. 



In this case, the concentration of occupied donors, 



n = N D - N A <^ N D . 



Thus all the electrons are able to find proper donors whose potential is lowered by the 
potential of charges neighborhood. Thus the Fermi level is situated well below —Eq, see 



Figure 3.4: Energy diagram of highly compensated semiconductor. Solid line shows the 
conduction band edge, occupied states are shaded. 

To understand the structure of the DOS tail consider a certain donor having at small 

1/3 

(r <C N D ) another donor which is positively charged . Its influence at r ^> a is just 
e = —e 2 /nr. It is assumed that the second donor is empty, the binding energy being ~ Eq. 
Thus a pair can contain only 1 electron. 

To find g(e) let us calculate the concentration of pairs having the energy in the interval 
(e, e + de). The probability to find another donor is AnN£,r 2 (dr /de) de, where r(e) = e 2 / ne. 
One has to multiply this quantity by No and divide by 2 (pairs!) to get 



Fig. M 



E 




g(E) 



9^) = l^N Dl ei) = ^ (47^/3) 1/3 . 

A 6 fa 



Now the Fermi level (at T = 0) is easily found from the conservation law 




— oo 



g(e) de = n . 




p 2 1 /3(i_;f)i/3 ' 
This is purely classical result. It is valid only at ry, ^> a. 
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Long-range fluctuations 

At 1 — K <C 1 it is a large and important effect because screening is weak. To obtain an 
estimates we can repeat the discussion for K <C 1 and replace the donor concentration 
by the total concentration, N t , of donors and acceptors. In this way we get for a typical 
potential energy of an electron, 

l(R) = ^(N t R^ . 

It also diverges at large R. 

How the screening takes place. The excess fluctuating density is 

AN = (N t R 3 ) 1/2 /R 3 - 

They can be neutralized at AN = n leading to the expression for the screening length, 

(7V tr 3)V2 2 2 / 3 iV7 1/3 

T) = — > r = 

r 3 r « (1-^)2/3- 

The random potential amplitude is 



nn 



1/3 



It increases with decreasing n. 

Here we disregarded the states with several extra electrons which are important very 
deep in the gap. 



3.3 Specific features of the two-dimensional case 

In high-mobility heterostructures the two-dimensional electron gas is separated from the 
doped region by a relatively wide undoped region, spacer, with thickness s. The fluctuations 
of the charged donors create a random potential i^(r) in the 2DEG plane. The reason 
to introduce a spaced is that the random potential becomes smooth and transport cross- 
section a small factor 2pps/h due to low-angle scattering. 

To calculate this potential it is practical to assume that the spatial distribution of 
charged donors in the doped plane is random and not correlated, 

(C( ri )> = , (C( ri )C(r 2 )> = C <J(n - r 2 ) . (3.6) 

For the field F b (r) we obtain 
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Here 

C(q) = j-^y Cij) e-^ -> (C( qi )C(q 2 )) = C(2n) 2 5( qi + q 2 ) . (3.8) 
As a result, we obtain 

(F 2 (r)> = 2ttC [°° ^V 2gs . (3.9) 

Here g m i n ~ 1/L where L is the linear dimension of the sample. Denoting 

e 2 , 

W = — V2^C (3.10) 

K 

we have 

y/{Ff) « W 0n(L/2s) . (3.11) 

For GaAs heterostructure with C = 1.0 x 10 11 cm -2 we have W = 9.2 meV. That means 
very large fluctuations of bare potential in a macroscopic sample. 

Electron screening removes this divergence. In the linear approximation it can be 
accounted for by replacing q by q + q s in the denominator of the integrand of Eq. (|3.9|) 



(F 2 (r)> = 2nC f°° -^-e" 2 ^ . (3.12) 

where 

e 2 dn , 

q s = 2n-—, 3.13 

K djJ, 

where n is the electron density while \x is the chemical potential of electrons. Following 
previous derivation we obtain 

W 

(^) = wr-^ for g sS »l. (3.14) 

We obtain that the linear screening dramatically reduces the potential since q s = 2/clb 
where is the Bohr radius. However, this estimate does not hold if the relative change 
in concentration is not small relative to the average concentration. Then the fluctuations 
in the long-range potential are of the order W. 

The main idea of nonlinear screening is that entire plane screens each Fourier harmonic 
of the impurity charge distribution independently. The excess average square density of 
impurities in a square of area R x R is VCR 2 /R 2 = \[CjR. If the electron density n is 
much greater, n 3> \/C / R, then the harmonic with the scale R will be screened linearly. 
However, the harmonics with R \[Cjn are not screened effectively, and 2DEG density 
becomes strongly inhomogeneous at the scale R. As a result, we arrive at the picture that 
harmonics in the impurity distribution with wavelengths R ^> R c = \fCjn are screened 
linearly, while those with R <C R c are screened very poorly. 

For large electron density, R c <^ s, and all significant harmonics in the bare potential 
are screened linearly. If the electron density is lowered, the nonlinear screening length 
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becomes larger, and eventually, for small enough electron density R c becomes of the same 
order as s. At small densities, the fluctuations include all harmonics between s and R c , 
and quantitative results can be obtained only numerically. We come to a conclusion that 
at some electron density a percolation between metallic droplets disappears, and a metal- 
to-insulator transition takes place. An of the critical electron density is 



n c (3 = 




(3 « 0.1. 



(3.15) 



s 



An extremely interesting problem appears in an external magnetic field. We will not discuss 
it here. 



CHAPTER 3. IMPURITY BAND FOR LIGHTLY DOPED SEMICONDUCTORS. 



Chapter 4 

DC hopping conductance in lightly 
doped semiconductors: General 
description 



4.1 Basic experimental facts 

At large temperatures semiconductors possess an intrinsic electrical conductivity due to 
thermal activation of carriers across the energy gap. In this case the electron, n, and hole, 
p, concentrations are equal and exponentially temperature dependent, be 

(27rkT^/m e m h f/ 2 Eg/2kT 
n ~P- 4yr 3^3 c • (4. ij 

Due to the activation energy E g /2 the concentration rapidly decreases with the tempera- 
ture. At low temperatures it becomes less than the concentration contributed by impuri- 
ties. In this region the conduction is called extrinsic since it depends on the impurities. 
A schematic plot of temperature dependence of the conductance is shown in Fig. |4.1| The 
region A corresponds to intrinsic conduction, while the regions B — D correspond to ex- 
trinsic conductance. If the impurities are shallow, then there exists the region B which is 
called the the saturation range. In this range, all the impurities are ionized and hence the 
carrier concentration is temperature independent. Consequently, the temperature depen- 
dence of resistivity is determined entirely by the carrier mobility. For example, a decrease 
in resistivity can be explained by a weaker phonon scattering. 

Further decrease in the temperature leads to a gradual freeze-out of the electrons which 
are recaptured by donors (for n-type material). 

Experimentally, the temperature dependence of the resistivity in the region C often 
shows only one activation energy, e±, 

p(T) = pe tl/kT , (4.2) 
where e\ is close to the ionization energy Eq of an isolated donor. 
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ABC d in 

Figure 4.1: Schematic temperature dependence of the resistivity of a lightly doped semicon- 
ductor. A - region of intrinsic conductance, B saturation region of impurity concentration, 
C - freeze-out region, D - hopping conduction region. 



At lower temperatures, hopping between the impurities without excursion to the con- 
duction band becomes most important (region D). The experimental results for p-type 
neutron-doped Ge is shown in Fig. |4.4 Each curve corresponds to a different acceptor 
(Ga) concentration. As a result, an interpolation formula for the regions C and D has the 
form 

p-\T) = pi 1 e~ ei/kT + P3 1 e~ e3/kT . (4.3) 

The 1st item corresponds to the band conductivity, and it is almost independent of the 
acceptor concentration. The second one corresponds to hopping conduction (it will clear 
later why we gave the subscript 3 for it). We see that that hopping conduction has 
a noticeable activation energy e%, though it is small compared to ex- It arises due to the 
dispersion of impurity levels. In hopping over donors, an electron emits and absorbs phonons 
that results in an exponential temperature dependence. 

When impurity concentration increases, it first somehow enhances that activation en- 
ergy. However, a further increase in the concentration increases an overlap between the 
wave functions of neighboring centers. At Na = 10 17 cm~ 3 the activation energy 63 van- 
ishes, and the conductance crosses over to a metallic one. 

A characteristic feature of hopping which is evident from Fig. [O] is an extremely strong 
dependence of on impurity concentration. It can be represented as 

Ps = P03 e /(JV -) (4.4) 

where p 3 an d f(N D ) are power-law functions of the impurity concentration. The expo- 
nential dependence is rather clear - it is just due to change of the wave-function overlap 
integrals. 

What is also important is that both p 3 and €3 are strong functions of the degree of 
compensation. The dependence £z(K) is non-monotonous, it decreases sharply at first, 
then has a minimum at K w 0.4, and then grows rapidly at K — > 1. 
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Resistivity of compensaied 
p-germanium for K - 0.4 J4,6], The 
curves correspond to the following 
acceptor concentrations (in cm" 5 ): 
(J) 7.5x10" 

(2) 1.4 x10 15 

(3) 1.5X10 1 * 

(4) 2.7x10 ,! 

(5) 3.6x10" 

(6) 4.9xl0 li 

(7) 7.2xl0 1i 

(8) 9xl0 ,s 

(9) 1.4 x10 16 

(10) 2.4xl0 lft 

(11) 3.5x10 ls 

(12) 7.3 xlO' 6 

(13) 1.0 xlO 17 

(14) 1.5 xlO 17 
(/5) 5.3 xlO 17 
{/(J) 1.35 xlO' 8 



0,2 0.4 0,6 0,8 

r' iK" 1 ] 



Figure 4.2: 



In addition to the band and hopping mechanisms of conduction , semiconductors with 
low compensation K < 0.2 display another activated mechanism which manifests itself 
near the Mott transition. The mechanism contributes one more term of the form 

p 2 l e~^ kT . 

We will not discuss this feature in detail. 



4.2 The Abrahams- Miller resistor network model 



4.2.1 Derivation 

To calculate the hopping conductance we shall use the so-called resistor network model 
introduced my Miller and Abrahams. First, starting from electron wave functions we shall 
calculate the hopping probability between the centers i and j with the absorption and 
emission of a phonon. Then we shall calculate the number of % — > j transition per unit 
time. In the absence of an electric field there is a detailed balance, and this number is 
exact equal to the number of reverse transitions. However, in the presence of the of a weak 
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electric field an imbalance appears, and a current appears proportional to the voltage drop 
between the centers. Consequently, we can introduce an effective resistance Rij of a given 
transition, and thus the whole problem is reduced to calculating of effective resistance of 
a network of random resistors. This way will be traced in detail below. 

Consider two donors having the coordinates and and sharing 1 electron. For 
simplicity, let us assume that the distance between the donors is large, so the overlap 
will be weak. In the absence of the overlap the wave function = ^(r — r^) satisfies the 
Schrodinder equation with the potential U = — e 2 /n\r — | . Within the framework of linear 
combinations of atomic orbitals (LCAO) method, the system is represented by symmetric 
and asymmetric combinations of atomic wave functions, 

= r («) 

v/2,/l±/dr***. 



The corresponding energies are obtained from the Hamiltonian 

,2 „2 



H = H - , - — r (4.6) 



K\r — Ti\ K\r — rj| 



are 

e 2 



E 2>1 = -E ±I ijt (4.7) 



where Eq is the energy is an isolated donor, while 7j, is the energy overlap integral is 



o2 



J. . = / *** 7 ~ r dr - I dr' I -. , dr . (4.8) 

J K\r — Tj\ J J k | r — r j | 

Now let us consider a simple case where the donor state is connected to a single extreme in 
the Brillouin zone center. Then we can split integral into a sum over crystal cells and use 
the fact that the envelope function is almost constant over a cell. The same can be done 
with potential energy which contains |r — r^j ^> a Q . The cell integrals |w ni o| 2 give unity. In 
this way we come to the expression 

J« = / FitfFJr) , e ' r dv - [ F^FAv') dr' f ^\ dr . (4.9) 



K\r — Tj\ J J n\r — Tj\ 



For a hydrogen-like function 



we get 



F(r) = —= e~ r/a 



ira 



3 




exp (-^j . (4.10) 

Note that there is a more accurate calculation which does not assume the overlap integral 
to be small. It differs from Eq. (|4.10|) by replacing 2/3 — > 2/e. The general feature of the 
energy overlap integral is that it can be expressed as 

1^ = 10 exp (-^) (4.11) 
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where 1$ is of the order of the effective Bohr energy and only weakly dependent on r^. 

In important fact is that the electron interacts also with the other charged impurities 
surrounding the center by some potential W(r). As a result, for the majority of donors, 



A{ > , where A{ = W{rj) - W(r< 



(4.12) 



Consequently, the potential W should be included into the Hamiltonian. The simplest 



situation case is when A] 3> L 



Then for the two lowest states we can show that 



\]>. + JL\u . 



A' 



A 



(4.13) 



Despite of the the fact that an admixture of the "foreign" state is small, it leads to the 
charge transfer on distance r^. 

Now let us consider phonon-induced transitions and assume that only a single branch 
of acoustic phonons is present. The transition probability i j with absorption of one 
phonon is 



2tt Vn 



\M^\ 2 5{hsq- A|)dq. 



(4.14) 



% (2tt) 3 

where Vo is the volume of the crystal, s is the speed of sound, q is the phonon wave vector, 
while 



M n = iD\ 



IhqN^ 
2V p 



%e ltl T % dr 



(4.15) 



is the matrix element of electron-phonon interaction. Here D is the deformation potential 
constant, iV q is the Planck function, p is the crystal density. Substituting the functions \&' 
and coming to the envelope functions we rewrite the previous expression in the form 



M n 



IhqN^ 
2V p 



\ 2" 



J if e iqr dr - J F f V q,r dr 



Ai 



Fji^e iqr dr 



(4.16) 



To estimate the integral let us evaluate two dimensionless parameters, 

qrij = Ajrij/hs and qa = Aja/hs. 



Since A? 



e^V 1 and R, 



N 



-y In 

D we see that qr^ is usually large (20-30) while qa is 
of the order 1. The characteristic length of variation of FiFj across the line connecting two 
centers is or the order Jr^a. Thus the last integral in Eq. (4.16) is strongly suppressed 
due to oscillatory factor, and it can be neglected. The first two integrals can be combined 
in the form 



[1 + q 2 a 2 /4) 2 
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The last equality is valid for hydrogen- like wave functions. Finally we obtain 



hqD 2 flijY iV q (l-cosq-r 
2pV s (1 + qW/AY 



Kl/ 



(4.17) 



Substituting this expression into (4.14) and using the fact that qr^ > 1 we get the final 
expression 



Hi = lij e 



p - 2r iil a 



N(A{ 



(4.18) 



where 



D 2 A 3 



N(A{ 



2e 2 V 



Trps 5 h 4 \3K,a J 
exp(A{/kT) - 1 



1 



[1 + {Ala/2hsfY 



(4.19) 



Let n = (0, 1) are the occupancies of ith donor which fluctuate in time. The transition 
% — > j is possible if = 1, nj = 0. Therefore the average number of electrons which make 
transitions is 

r - = (7^(1 - n^)) . (4.20) 

The average should be taken with respect to time. Here we make a fundamental assumption 
that the occupancies can be replaced by their average numbers, (n^) = fi, while the energies 
can be replaced by self-consistent energies 



acc 



don „2 



— i 

l K \ r i — r l\ 



E 



e 2 (l-h) 

lA T % ~ r k\ 



In this approximation we can put A{ = Ej — Ej = to obtain, 

\'u ~u< 2r '"Mr J -:;)f;(\-f i ). 
For the reverse process we get 

I7i = 7° e- 2 ^>' [iVfe - £l ) + 1] fjil-fi) 
The current between the centers j and i can be written as 



J i 



-e{Tij - r 



In the absence of an electric field 

fi ~ fi 



cxp 



so - V 
kT 



+ 1 



(4.21) 



(4.22) 



(4.23) 



(4.24) 



(4.25) 



where Eq is the average energy at the site at E = 0, while factor 1/2 enters because there 
are 2 possible spin states at the site i. In fact the energies Eq and the functions ff are 
made self-consistent by Eq. Q4.21 ). 



4.2. THE ABRAHAMS-MILLER RESISTOR NETWORK MODEL 



45 



It can be easily seen that a detailed balance is present, and there is no current. The 
electric filed, first, redistributes the electrons between the donors, jf — > ff + 8fi. We 
describe these changes by small additions 8 fa, 



1 fe - 5 fa 
2 GXP W 



-\ -i 



+ 1 



(4.26) 



Secondly, the field affects the donor level energies, 



g 2 don 

Si = £° + 8ei , 8ei = eE • rH ^ 



Sf k 



r k \ 



(4.27) 



Would the circuit be broken, a new equilibrium with 5 fa + 8ei = will result. Howveer in 
the case of conductance, 8 fa + 8ei ^ 0. Assuming all the correction to be small {linear in 
E effects) we can expand all the energies up to the first order in 8ei and 8 fa. After such 
manupulation we obtain 



J, 



er°. 

-j^r [8 fa + Set - (8 fa + Sef)] 



The latter equation can be expressed as the Ohm's law 

J ii = RZ 1 {U i -U J ) 

where 



Rij 



kT/e*T% , Y% = 7 ° e" 2 ^ N(e 3 - e t ) - /°) , 

„2 don 



-et/j = <5/ii + 8ei = eE • H ^ 



- r fc | 



(4.28) 

(4.29) 

(4.30) 
(4.31) 



The quantity — eUi can be regarded as a local value of the electrochemical potential on 
the donor % counted from the chemical potential fa Thus Ui — Uj is the voltage drop on 
the transition % — > j while R^ is its resistance. We end at the problem of calculating of 
effective resistance of a random resistor network keeping the potentials at the electrodes 
constant. 



4.2.2 Conventional treatment 

Let us analyze the expressions for the resistances R^ and assume the temperatures to be 
low, 

fcT«| £ °- S °|, 

Then we can extract the functions strongly dependent on the distance and on the 
energies, 

r°- ~ lij e~ 2rij/a e~^ /kT (4.32) 
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where 

e ij = g (\ £i ~ £ i \ + \ £i ~ A*l + kj - ■ ( 4 - 33 ) 
Now it is convenient to introduce the quantity 

& = 2 ? + I < 4 - 34 > 

to express random resistors as 

i^ = 4exp(^), R% = kT/e 2 ^ . (4.35) 

The main feature of the problem is an exponentially wide range of resistances. It is clear 

that doth averaging of conductances and resistances will lead to wrong answers. A correct 
method of averaging is suggested by the mathematical theory of percolation. 



Chapter 5 
Percolation theory 



The term "percolation" has been introduced in 1957 by Broadbent and Hammersley in con- 
nection with a new class of mathematical problems. We shall start with general description 
and then will map them to the hopping conductance. 

5.1 Lattice problems 
5.1.1 "Liquid flow" definition 

Imagine an infinite lattice with bonds between the adjacent sites which permit the flow in 
both directions, and that each wet site will instantly wet the adjacent sites. The various 
problems follow from introduction random elements into the theory. 

Bond problem 

Assume that each bond in the lattice can be other blocked, or unblocked. Let the proba- 
bility of any bond to be unblocked is "x" . The distribution of unblocked and blocked bonds 
remains constant in time. 

If a random site is wet, then there are 2 possibilities. The initial site will wet either 
a finite or an infinite number of other sites. The key quantity is the probability of a 
random initial site to wet an infinite number of other sites. In the infinite lattice this 
probability depends only on x. Let us denote it as P b (x). The bond problem is illustrated 
in Fig. |5.1| (left panel), while the function P b (x) is plotted for different lattices in the right 
panel. Curves 1,2 and 4 correspond to 3D lattices, the rest to 2D ones. When x is small 
P b (x) = since the blocked bonds prevent the liquid to spread far from the initial site. 
On the other hand, as x approaches 1, so does P b (x). A very important quantity is the 
percolation threshold, x x , defined as an upper limit for the x-values for which P b (x) = 0. 
When x-x c < 1, one has 

P b (x) ~ (x - x c f (5.1) 
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Figure 5.1: 



where b is a critical exponent. The situation resembles type II phase transitions. Note that 
all the conclusions are valid for an infinite lattice. 

One can also define the probability to wet iV other sites, P^(x). P^{x) differs from 
at < x < 1, but at N ^> 1 it is small for x < x c . We can define 

P\x) = lira P b N {x) . 

An example of the bond problem - an orchard which is planned in the form of a square 
lattice. It is known that a diseased tree contaminates another tree away at the distance 
r with the probability f(r). It is required to find a minimum lattice period h m i Q to avoid 
epidemy. The solution is obvious, 

f{h min ) = x c . (5.2) 



Site problem 

For the site problem, sites can be either blocked, or unlocked, while the bond are assumed 
to be ideal. Blocked sites permit no flow in any direction, they cannot be wet. Let x be 
the fraction of unblocked sites while P s (x) be a probability for a random site to wet infinite 
number of other sites. 

The probability function P s (x) is plotted in the left panel of Fig. |5.2| , while a realization 
of the site problem in shown in the right panel. One can see that there is a percolation 
threshold which depends in the lattice type. 



5.1.2 Cluster statistics 

Both the site and both problems are usually discussed in terms of cluster statistics, rather 
in terms of liquid flow. Let us apply this approach to the site problem. Imagine that a 
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Infinite cluster density 
P {s} (x) for three cubic lattices [5.4J: 
(/) simple cubic, (2) body centered 
cubic, (J) face centered cubic 




Site problem on a square lattice 



Figure 5.2: 

fraction x of the sites is painted black, while the rest sites are painted white. Any two 
adjacent sites are connected if both are black. A chain of such connected clusters becomes 
a cluster. For example in the right panel of Fig. there are 3 clusters - one consists of 5 
black sites while two consist of 3 black sites. 

When x is small, then there are only small clusters. As we approach the percolation 
threshold, several clusters may come together, and the mean cluster size increases. At 
x —>■ x c an infinite black cluster is born. It resembles a random network which permeates 
the entire lattice, while smaller since finite clusters exist in the "pores", see Fig. |5.3| . 
Actually, P s (x) has a meaning of the density of the infinite cluster. Similarly, a proper 




Schematic representation 
of infinite cluster and several 
large (critical} finite clusters for 
< x- x c 1 . The lines are actually 
chains of black sites, but ihc resolu- 
tion chosen is so low thai details of 
the order of tbe original lattice period 
are not noticeable, L(x) h the corre- 
lation length. The dashed squares 
have sides !< L(x) and t>L(x) 



Figure 5.3: 
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definition for the bond problem can be formulated. 

It is assume that no more than one infinite cluster can exist in a lattice. Than brings 
us to an important application of the site problem - a crystalline solution of ferromagnetic 
atoms (with a fraction "x") in a non-magnetic host. Assuming nearest-neighbor ferromag- 
netic interaction we arrive at the ferromagnetic phase transition at zero temperature as a 
percolation threshold for the site problem. Consequently, 

• the transition takes place at x — x c ; 

• the saturation magnetization at T = for x > x c can be expressed through the 
density of the infinite cluster, 

M\ T=0 = fi P s (x) (5.3) 
where fio is the atomic magnetic moment. 

Another important application is the aforementioned Anderson transition. The appearance 
of the extended states is similar to formation of an infinite black cluster. 

Now let us introduce relevant quantities for the cluster statistics. Let n s be the number 
of clusters containing s sites, per site of an infinite lattice. Then the fraction of sites 
belonging to an infinite cluster is 

P{x) = x — sn s . (5.4) 

s 

The second important quantity is the mean cluster size, 

S{x) = — 

According to numerical estimates, S(x) becomes infinite at x c , 

S(x) ~ \x — x c |~ 7 , (\x — x c \ <C 1) , (5.6) 

where 7 is a critical exponent. Note that Eq. calculates a weighted average which is 
divergent. A simple average would be regular because 

lim sn s = x c . 

S 

The majority of sites remain in clusters with s ~ 1 even at x = x c . 

According to numerical estimates, n s near the percolation threshold behaves as follows. 
There is a critical number s c of the sites in the cluster, which grows at x — > x c + and at 

s c (x) ~ \x - x c \~ A , A>0. (5.7) 

If s s c and s increases the quantity n s decreases according to some power law, while at 
s ^> s c it decreases exponentially. Consequently, the numerator of Eq. ( |5.5|) is determined 



(5.5) 



5.1. LATTICE PROBLEMS 



51 



by the clusters with s ~ s c , which we will refer to as critical clusters. The divergence of 
S(x) means that the number of sites in critical cluster increases. 

The third important quantity in the percolation theory is the correlation function. Let 
us define a function g(ri,rj) as equal to 1 if both sites are black and belong to the same 
cluster, and otherwise. Now we can introduce the pair correlation function as 

G(r, x) = G{vi - rj, x) = (g(r h r,-)) . (5.8) 

Naturally, G(r,x) does to zero at r — ► oo because n s decreases with s. At the distances 
smaller that the average size of the critical cluster which we will denote as L(x), the 
correlation is governed by the clusters with s <C s c (x). Consequently, G(r,x) decreases as 
a power law of r. If r 3> L(x), the exponentially rare clusters with s ^> s c (x) are dominant, 
and G(r, x) decreases exponentially. The length to the critical clusters, L(x), is called the 
correlation length. Since the number of sites in the critical clusters grow at x — > x c ± 0, 
the correlation length must grow if one approaches the percolation threshold, 

L(x) ~ \x — x c \~ u , (\x — x c \ <C 1) . (5.9) 

Here v is the critical exponent of the correlation length. It is usually thought that the 
critical exponents 7, A and v are the same whether x > x c or x < x c , and that there exists 
a sort of cluster-pore duality (see Fig |5.3|) . Thus, at x < x c the quantity L{x) is a typical 
radius of clusters, while at x > x c it is a typical radius of the pores. Near x c the correlation 
length introduces a typical distance ( "period" ) between the nodes of the network. 

5.1.3 Percolation through a finite lattice 

Consider a finite lattice containing I sites per side. Let all the sites be initially white. We 
will randomly paint some sites black, slowly increasing the fraction x. At some x = x^ the 
black cluster will connect opposite sides. We will consider percolation threshold 

for a finite lattice. If we permit the procedure, the arrangement of the black sites, as well 
as x c i will be different. Therefore, x^ is the random quantity, and we can define the true 
percolation threshold as 

x c = lim (x ci ) . (5.10) 
£^00 

Another variation of the same theme is to define x as the probability for a site to be black 
rather that as the fraction of the black sites. In a finite systems, these two definitions lead 
to different valuers of x c e and (x c i) due to fluctuations. We can say that the 1st definition 
corresponds to a micro-canonical ensemble while the 2nd one - to the canonical ensemble. 

Figure f[Q shows the distribution function of Xd in a canonical ensemble for the site 
problem on a square lattice. This figure was obtained by Monte Carlo simulations. The 
distributions are nearly Gaussian in shape. The larger the side the lower the dispersion 
We = ({x c e — (x c i)) 2 ). It was shown that We obeys the power law 



W e = B£- 1/v , 



(5.11) 
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The distribution functions of values of ihe percolation threshold for the site problem on 
/x / square lattices [5.13]. The / values in units of the lattice constant are: (/) 8, {5) 32, (J) 128 

Figure 5.4: 

while the average threshold obeys the law 

( Xct ) = x c + Ar 1/X . (5.12) 

In three dimensions, v m 0.9 and A ~ 1. It can be shown that the exponent v is the same 
as for the correlation length. 

As an example of percolation between the opposite faces for the 3D bond problem let 
us consider electrical conductivity of a large cube in which the unblocked bonds are repre- 
sented by standard resistors. The blocked bonds are assumed to have infinite resistance, 
Fig. |5.5| . If the cube is large, then the conductivity differs from zero only at x > x c (b), 
while at < x — x c <C 1 one has 

a{x) ~ (x - x e f , (5.13) 

where t is yet another critical exponent. It is clear that to find a it is not enough to know 
P b (x) since also a topological structure of the infinite cluster is important. 

5.1.4 Summary 

We have presented 3 definitions of the percolation problem: liquid flow, cluster statistics, 
and percolation through a finite lattice. It can be shown that critical phenomena are the 
same in all three formulations. The typical questions to be answered: 



5.2. CONTINUOUS PROBLEMS 



53 



n 



Figure 5.5: A random resistor network between two electrodes. 



1. What is the percolation threshold x c ? 

2. What are the critical exponents /3, 7 , v and tl 

3. How do dimensionality and lattice geometry affect the above mentioned quantities? 

4. What is the topology of finite and infinite clusters? 

The values of the percolation thresholds are summarized in Table ^TT| In the Table ^2] the 
most important critical exponents are presented. The majority of values in the tables were 
found by numerical simulations, however some results are obtained analytically. In rare 
cases, experimental methods were used. 

The modern point of view is that the critical exponents depend only on the dimension- 
ality of the problem, and in this way they are universal. 



5.2 Continuous problems 



Suppose a random continuous function V^(r) is defined in the entire space by its correlations, 
such that (V(r)) = 0. Then let us paint black the regions where V^(r) < V and leave other 
regions white. As V changes from —00 to 00 the volume of black region changes form zero 
to all space, see Fig. (STfil . As V increases, the black regions grow, and at some V c they 
merge into an "ocean" The quantity V c is called the percolation threshold. Such a "liquid- 
flow" definition can be easily reformulated into the "cluster statistics", or "finite- volume 
percolation" ones. 

One can easily map the continuous problem to the lattice one. Assume a lattice with 
so fine grid that the potential does not change. Then fix V and consider all the sites in 
the black regions to be black while the other sites - to be white. The the fraction of the 
black sites is 

x = 0(V) = f F(V')dV (5.14) 
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Percolation threshold for 


various lattices. Exaci values are underlined. All oilier 


values were obtained by numerical methods 
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Table 5.2: From the book [1] 
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where 6(V) is the space fraction of the black sites, while F(V) is the distribution function 
of the quantity V. The percolation threshold then can be defined by equating the quantity 

9 C = [ Vc F(V)dV (5.15) 

to the percolation threshold x c of the site problem. Physically it means that the distri- 
butions of black and white sites are now correlated according to the function V(r). This 
proposition is very important for practical calculation of 8 C . 

Along with the above mentioned "percolation through black" , one can define a "perco- 
lation through white" with some threshold V' c . In the two-dimensional case V c — V' c — 0, 
and 6 C =0.5. The situation in 3 dimensions is more complicated. Only an inequality 
V c < V'c does exist, so the following 3 situations are possible: 

1. V <V C and there is percolation through white, but not through black; 

2. V c < V < V' c and there is percolation through both black and white; 

3. V > V' c and there is percolation through black, but not through white. 
It the potential is statistically symmetric, then 

V c < , V' c > , and 9 C < 0.5 . (5.16) 

The quantity 8 C has been extensively studied for various distribution functions of the 
random potential. A typical value of 9 C for Gaussian potentials is # c =0.17. 



5.3 Random site problems 

The most important problems for hopping conductance are the so-called random-site prob- 
lems. Random sites are the sites chaotic distributed in space, the average number iV of 
sites per unit volume being known. If one introduced some function £ii( r ii)j the two sites 
are considered bonded if, for some number £, 

C,< C (5.17) 



Eq. (|5.17|) is referred to as the bonding criterion. If two sites are bonded either directly, 



or via some other sites, then they form a cluster. It is required to find the percolation 
threshold £ c , i. e. the lowest bound of the parameter £ which still permits an infinite 
cluster. 

The simplest example of the bonding criterion is 

rij<r, (5.18) 



which is met if the state j is located within the sphere of the radius r centered at the site 
i. Consequently, the problem has a very simple geometric interpretation: spheres of radius 
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r are constructed around all the sites, and we must find the lowest value r = r c which 
allows endless chain of sites in which each cite lies inside the sphere about the preceding 
site, see Fig. |5.7| We will term the quantity r c the percolation radius. Occasionally it is 
more convenient to construct spheres with radius r/2 and calculate the value that allows 
for an infinite cluster of overlapping spheres. 

The percolation threshold r c depends solely on the concentration N and is must be 
proportional to iV -1 / 3 . In most cases people consider the threshold value 

B c = (4/3)vriVr 3 (5.19) 

of the dimensionless parameter B = (4/3)irNr 3 which corresponds to an average number 
of bonds per site. This problem is often called the sphere problem, or the disk problem in 
2D case (where B c = irNr^). 

A physical implementation of this problem is a ferromagnetism of a dilute crystalline 
mixture of magnetic and non-magnetic atoms, where the exchange forced between magnetic 
atoms are modeled by a finite-radius interaction with the interaction radius r tnt . If r int is 
much greater than the lattice constant so that many sites fall into the sphere of interaction, 
then we arrive at the sphere problem. However, the problem is a reversed one: knowing 
r int we have to find N. If the number of sites that fall into the interaction sphere is Z that 
for the lattices of any dimensionality we obviously have 

lim Zx c = B c . (5.20) 

Z^oo 

The critical indices for the random site problems can be introduced in the same way as 
for the lattice model. Introducing n s as the number of clusters of size s in a given volume 
we can rewrite scaling properties in a similar way in terms of £ — £ c . 

In some important cases the bonding criterion is fulfilled on a complex surface rather 
than on a sphere. Then it should be written in the form 

C ; -ofr,,)< £. (5.21) 

Here £ > 0, and <j>(r) is a monotonously increasing homogeneous function of x^, yij, and 
Zij. The equation <p{rij) = £ define a surface which replace the sphere in the constructions 
discussed above. In particular, we can define B c as NV^ C . It can be shown that the values 
B c are the same for all the surfaces which can be obtained from each other via a linear 
coordinate transformation, 

X' k = ^ A k £X£ 

i 

since such transformation does not change the topology of the system. Thus the threshold 
values B c for ellipsoids and spheres are the same. 
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5.4 Electric conductivity of random networks of con- 
ducting bonds. Infinite cluster topology 

As we have mentioned for the case of lattice models, the electric conductivity of a lattice 
of conducting and non-conducting bonds behaves as 

a(x) ~ (x - x c y (x>x c ). (5.22) 

The aim of this section is to discuss the critical exponent t. There are few ways to study 
this problem: an experimental (by punching holes in a graphite paper, cutting bonds in 
a metallic grid, etc), computer simulations, etc. It appears that the results both for site 
and bond problems are closed in value. It seems logical that these values will retain also 
for the continuum version. What we want to do now is to relate the conductivity to other 
properties of the infinite cluster, namely, with its density P(x) and the "period" L(x). It 
is clear that these quantities do not contain the full necessary information about cr(x), and 
we need some assumptions to relate the critical exponents. 



Dead ends 

Consider a large cube containing I units per side. If the infinite cluster resembles a cubic 
lattice then 1/3 of the bonds form chains that connect the opposite edges of the cube, the 
number of bonds in these chains being « (l/3)P(x)£ 3 . Each chain has £ sites, therefore 
the number of parallel chains is ~ (l/3)P(x)£ 2 . Since the resistance of each bond is oc £, 
the resistance R of the whole cube is oc £[P(x)£ 2 ] -1 . Consequently, 

a{x) ~ {MY 1 ~ P{x) ~ (x -x c f. (5.23) 

We arrive at the conclusion that t = (3. However, we know that t > (3, see Table ^T2|- That 
leads to a conclusion that almost all the "mass" of the infinite cluster is concentrated in its 



dead ends, sii Fig. |5.8| . These chains contribute to P(x), but not to a(x). The insufficiency 



of the infinite cluster can be also due to redundunt 



The nodes and links model 

This is the most used model (Skal-Shklovskii-De Gennes, SSDG) assumes that there is a 
backbone network with w mean size L ~ (x — x c )~ u , and that this back bone cartries dead 
ends. The part of the backbone is called the link. It is assumed that at least half of the 
links are not doubled. Each link can be twisted, so its length L is larger than L. It this 
way we can write at x — x c <C 1 

L ~ (x — x c )~^ , where £ > v . (5-24) 

A schematic representation of the SSDG model is shown in Fig. |5.9| . To proceed let us first 
show that £ = 1. Let x > x c and cut each bond with the probability (x — x c )/x. As a 
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result, the fraction of unbroken bonds will become x[l — (x — x c )/x] = x c . Consequently, 
we arrive at the percolation threshold for the original lattice. Now notice that to cut a link 
we must cut one of its constituent bonds (we have assume that the links are not doubled!). 
Thus, the probability for a link to be broken is proportional to L and equal L(x — x c )/x. 
We can assume the percolation threshold for a backbone to be equal to some y c . Equating 
L(x — x c )/x = y c we find that ( = —1. We note that in 2D case v > 1, and it means that 
SSDG model does not hold, and that duplication of links is crucially important. However, 
for d = 3 z/3 ~ 0.9, and L ~ L 11 . Note that the twistness appears very weak comparing to 
the cases of the random walk, L ~ L 2 and of the self-avoiding random walk, L ~ L 5 / 3 . 

Now let us calculate a(x) If resistance of each bond is R then the resistance of a cube 
with side L is Rq(L/ 1$) where £q is the bond length. At the scales larger than L the infinite 
cluster is homogeneous. Therefore a(x) can be calculated as an effective conductivity of a 
cube with the side of the order of the correlation length. We get 

a(x) = R- l L 2 - d = Rq 1 (L / L 2 ~ d = R^H 2 ^ d {x - x c ) l+v{d - 2 ^ . (5.25) 

As a result, 

t = l + u(d-2). (5.26) 

For d = 3 we get t 3 = 1 + u 3 . Substituting z/ 3 = 0.9 we obtain t 3 = 1.9 while the numerical 
value is 1.7. We conclude that SSDG model is not very accurate. However it is extensively 
used because it is simple and qualitatively correct. It has been extensively used in the 
theory of dilute ferromagnets, granular superconductors, conductivity and Hall effect in 
disordered systems. The model has been also extended to take in some way into account 
the bond doubling. It makes the model more accurate. 

The most important extension (Pike and Stanley) is to discriminate between red and 
blue backbone bonds. Breaking a red bond breaks the link itself. It was further suggested 
that blue bonds from blobs connected by singly connected chains. The number of red bond 
between the nodes has also be scaled as L' ~ (x — x^ 1 . A finite size of the blobs removes 
the limitation of the model for the cases when v > 1, in particular in 2D case. It seems 
that the combined nodes-links-blobs model describes properly the structure of the infinite 
cluster. 

Scaling hypothesis and calculation of a(x) 

Another way to estimate the exponent t is to suggest that the large-scale structure of 
the network remains self-similar as the system approaches the percolation threshold. The 
second assumption to be maid is that the resistance is determined only by a large-scale 
structure. The two assumptions allow one to relate t and v in a simple way. 

Let us denote r = x — x c and examine the systems at two values of r, such that r 2 < t±. 
Since L ~ M^, the dimensions of the system will increase by the factor L(r 2 )/L(r 1 ) = 
(r 1 /r 2 ) t '. We assume that we also scale the distance between the contacts and the potential 
difference by the same factor. We conclude that the current distribution in the bond will 
not be affected. Now consider an area element that is perpendicular to the average direction 
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of the current. The number of channel piercing the area will decrease as L(ti)/L(t2) = 
{J2/t~i) 2u ■ Thus the current and the conductivity will decrease by the factor (ri/r 2 ) 2y . As 
a result, in three dimensions t = 2u, while in two dimensions t = v. In this way we obtain 
the values £3 = 1.80 and t<i = 1.33 close to the experimental ones. The scaling approach 
has appeared extremely fruitful during last years and revealed fractal properties of the 
infinite cluster. 



5.5 Percolation theory and electric conductivity of 
strongly inhomogeneous media 

We have previously shown that calculation of hopping conductivity can be reduced to 
conductivity of a random network with an exponentially wide range of resistances. The 
similar approach is valid also for a continuum. 

Suppose that we have a medium in which local conductivity fluctuates as 

d(r) = ffe^ r) , and (5.27) 
(K(r)-<0] 2 >»1- (5-28) 

For example, the source of fluctuation can be a band bending due to a random large 
potential, and £(r) = e(r)/kT. Here e(r) is the distance between the Fermi level and the 
bottom of the conduction band. A similar problem can be also formulated on a resistor 
lattice with R = R e^ where the random quantity £' is randomly distributed between — £ 
and £ > 1- 

To illustrate the method of conductivity calculation let us consider a cubic lattice. Let 
us first derive the exponential factor in the effective conductivity. Then we will discuss 
also a pre-exponential factor, and finally compare the result with the SSDG model. 

Consider a simple cubic lattice and select a value £ in the interval (— £o,£o)- Then let 
us replace all resistances with £' > £ by infinite resistances thus breaking the circuits. Let 
us call the so-obtained conductivity as cr(£). The chosen £ determines the probability of a 
random resistance not being broken: 

x(0 = f F{g) dg , (5.29) 

where the distribution function F is defined as 

f( « - { T°' II > !: • (5 - 3o) 

From Eqs. (^29|) , (|O0|) we find 

x(0 = ^- (5.31) 

Now let us increase the value of £ from — £ till some value £ c corresponding to the per- 
colation threshold for the bond problem, x c (b), given by Eq. (|5.31|) . If we that increase £ 
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from £ c to £ c + 1, the presence of infinite cluster will cause a rapid increase in conductivity. 
This happens because the correlation length rapidly diminishes as 

HO ~ io[x(0 - x c (Or ~ toCoti - Q~ v • (5.32) 

Therefore, the number of parallel conducting chains in the infinite cluster network will also 
rapidly increase. On the other hand, the change in the individual resistance remains small 
if £ is changed less than 1. Than implies a rapid power-law increase in a, 

f(0 ~ (6 - Q b , where b > . (5.33) 

Let the infinite cluster that is formed when £ and £ c differ by about one be called the critical 
subnetwork. Its resistance is determined by the highest resistances because, by definition 
these resistances cannot be shunted by lower ones since that would allow for percolation 
a t 6 < 6c- Consequently, for s critical subnetwork we have 

a(£ c + 1) ~ tree"*" . (5.34) 

Further increase of £ will not raise a significantly, even though the density of infinite cluster 
will rise. Indeed, if we continue to increase £ we introduce exponentially large resistances. 
Thus, new chains will be shunted by the critical subnetwork. The expression 



a = a, 



<c (5.35) 



is the foundation for the hopping conduction theory. According to the Table |5.2| , x c (b) 
0.25 for the simple cubic lattice. From ( |5.31|) we find £ c = — £o/2, and 



a = 4° /2 (5.36) 

This formula is valid only if £ ^> 1, when the concept of the critical subnetwork is valid. 
Numerical tests show that the model is practical at £o > 9. 

Let us now discuss the pre-exponential factor ctq. From dimensional arguments, 

v = c(to)Wor 1 (5.37) 

where c(£o) is a dimensionless factor which can contain only powers of £o- It can be shown 
that this power equals to — v where v is the correlation length exponent. 

To obtain the power, let us suggest another derivation of the conductivity. Consider 
a cube with side i cut from the infinite lattice. Let us break all the resistances and then 
connect them back in the ascending order until we reach the percolation between the 
opposite sides. Let the last (greatest) resistance be 6^. Therefore, 

R = R & -> -> o = (iV) _1 e~^ • 

The quantity is actually random, however for a large system its fluctuations are small. 
As we have seen, the fluctuations behave as 
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Now let us increase the cube size t. At some size the conductivity becomes I- independent. 
We can tell that it happens if the mean square fluctuation approaches unity, or at 

Lq is just the infinite cluster correlation length. For such a cube the difference between £^ 
and £ c is negligible. Thus, 

A same expression can be easily derived from the node-link model of the infinite cluster. 
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Percolation in a two- dimensional potential 
1 5. 59]. The outlines correspond to surface 
V{i) = const. "Black" regions fl^(r) < V) are shaded. 
The three maps correspond to three different values of 
V{V X <V 2 < V 3 ) 



Figure 5.6: 
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Figure 5.9: Schematic representation of the SSDG model. 



Chapter 6 

Hopping conductance 



6.1 Dependence on impurity concentrations 

In this section, we shall apply the percolation theory to study the dependence of the ps 
contribution on the impurity concentration. 

6.2 The contribution p% as a function of impurity con- 
centrations 

Let us concentrate on the simplest case of isotropic impurity wave function. In fact, we 
have to calculate the resistance of the Abrahams-Miller network of random resistors, 

i^ = i2° e««, (6.1) 

where 

& = 2r ij /a + e ij /kT, (6.2) 

R% = kT/e 2 ^, (6.3) 

€ij = (1/2) (\€i - €j\ + \€i - (Ji\ + \€j - (Jl\) . (6.4) 

To calculate the resistance let us apply the percolation approach. First suppose that only 
the resistances with £y < £ are switched on. Let us gradually increase £ until we reach 
the condition of percolation over conducing resistances. That obviously occurs at £ = £ c . 
Increasing £ from £ c to £ c + 1 results in formation of the critical subnetwork of resistances. 
This subnetwork shunts all the resistors for which £ > £ c +l. Consequently, the exponential 
part of the hopping conductivity to the parameters of the problem is given by 

p = Po e 5c . (6.5) 

This is the starting point of the theory of hopping conductivity. 
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Let us discuss the case of rather high temperatures when we can neglect the item 
6ij/kT in Eq. (|6.4j ). The bonding criterion has now a very simple form, < a£/2, i. 
e. the connectivity of a pair depends only on their separation. In this way we return to 
the problems of spheres discussed in the previous chapter, see Fig. |5T7| . Therefore £ c can 
be expressed in terms of the percolation radius r c , i. e. the percolation threshold for an 
auxiliary problem with the bonding criterion < r as 

& = 2r c /a . (6.6) 

According to numerical calculations, 

B c = (An/3)Nr 3 c =2.7±0.l, r c = (0.865 ± 0.015)iV" 1/3 . (6.7) 



Here N is the concentration of the majority dopant. Substituting (|6.7|) into (|6.6|) , we find 
from (|6.5|) 

p 3 = p 03 exp[a/{N 1/3 a)] , a = 1.73 ±0.03. (6.8) 

This result agrees with the experiment. The problem to calculate the activation energy is 
more delicate and the answer depend strongly on the degree of compensation. 

The case of anisotropic wave function is different just geometrically because the bonding 
condition becomes anisotropic. 

A very interesting changes occur in an external magnetic field. The presence of mag- 
netic field squeezes the electron wave functions, and the overlap integrals between the 
components of the pair decease. Furthermore, the wave functions are squeezed mostly in 
the direction perpendicular H. As a result the resistivity exponentially increases in a strong 
magnetic field, the resistivity becoming anisotropic. Giant positive magnetoresistance is 
the hallmark of hopping conductance. 



6.3 Activation energy 

6.3.1 Low degree of compensation 

As we have seen,the resistivity of a lightly doped semiconductor in the hopping regime has 
the form 

p = p 3 e^ kT . (6.9) 

Since we were interested only in p 3 the energy term e^/kT has been ignored in the ex- 
pression for t^ij. To take it into account let write £y in the form ^ = £?■ + where 
£?• = 2r ij /a while A£ y = e^/kT. 

Let us split the percolation threshold into two parts, 

e c = e c o +A£ c (e.io) 

where £° is the known percolation threshold for the problem with bonding criterion 

$ ^ * ( 6 - n ) 
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while A£ c is the correction which we are going to find. 

In the case K <C 1 let us first ignore the long-range potential. Then the most donor 
levels are close to the unperturbed value (which according to our notations corresponds to 
e = 0) and well separated from the value of chemical potential. Then for overwhelming 
majority of donor pairs we can neglect q and €j in comparison with p. As a result, 

A^^p/kT, (6.12) 

and 

£ c = e c +v/kT. (6.13) 

As a result, 

p = Po3 e ^ kT , e 3 = 0.99e 2 N^/K . (6.14) 

Since the long-range potential is small comparing to p it produces only a small correction 
to the activation energy, 

e 3 = ^(i - 0.29fT 1/4 ) . (6.15) 
6.3.2 High degree of compensation 

Neglecting the long-range potential we can repeat all the estimates for the case K <C 1 
and set = \p\ . In this way we arrive at the conclusion that 

e 3 = M = eD2- 1/3 (l-^)- 1/3 . (6.16) 

However, in most cases the long-range potential appears important. Since the long-range 
potential produces smooth distributions, in the vicinity of a point r the donor energy is 
mainly determined by the energy 

V(t) = &j>(r) . (6.17) 

The neighboring donors have almost the same energies (smooth potential!!). So, since 
^( r ) > P, w e get 

p(r) = p 3 e ^W-M]/w . (6.18) 

Now we arrive at the problem to calculate effective resistance of the system with an ex- 
ponential spread of resistances, and we can apply the percolation theory for the second 
time. 

Let us paint black the regions where V(r) / kT < £ leaving the remainder white. Suppose 
that only black regions are conducting while the white ones are insulating. Now begin 
increasing £, thus including the regions with higher resistivity. At some £ = £ c = V c /kT 
the black regions merge into an infinite cluster, and further increase in £ does not effectively 
change the resistivity. Therefore the resistivity is determined by the critical subnetwork 
and given by the expression 



P = P3 e^ +v ^ kT , -> es = V-p . 



(6.19) 
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The quantity V c can be found by numerical calculations. It is of the order of a typical 
fluctuation in the potential energy, 

e 3 = C l7 (r s ) = C 1 e D (l-K)- 1 /3 ) Cl ~ 1 . (6.20) 

Form the above consideration we find the law 

ei = e 3 + E (6.21) 

which relates the activation energy to the conduction band to those in the impurity band. 
This relation has been confirmed by many experiments. 



6.4 Variable range hopping 

At low temperatures hopping occurs between the states whose energies are concentrated 
near the Fermi level. Because of the factor exp(ejj) / kT in the expression for the resistance 
Rij, only the states within the narrow band near the Fermi level are important. The 
quantity are given by Eq. (|4.32|) . At sufficiently low temperature the density of stated 
in this band can be regarded as constant, g(e) = g(f-i), see Fig. |B~T| . Similarly, we assume 

a £ 
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-e- 



-e- 



-e- -e- 
-e- 
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Figure 6.1: Construction of the band containing the states within the layer /i — eo, p — e . 



that the localization radius is also energy-independent. 
Let us consider the contribution of a small energy band 



(6.22) 



which is symmetric about the Fermi level, see Fig. |6.1| . The concentration of states in this 
band is 

N(e )=2g(p>)e . (6.23) 

The typical resistance can be obtained by replacing the distance by the average one, 
[A^eo)] -1 / 3 , and — > eo- Then we have 



p = po exp 



1 



e 



[^(eo)] 1 /^ kT 



p exp 



e 



[g(p)e ]^ kT, 



(6.24) 



6.5. COULOMB GAP IN THE DENSITY OF STATES 



69 



Here we ignore numerical factors. Let us analyze the dependence p(eo). It has a sharp 
minimum at 

{kTf^ 



£0 " C0(T) " [ 9 (pW}^ ■ 
Substituting this value into Eq. Q6.24J) we get the famous Mott's law. 

P 



P(T) = Po e^ 1/4 , T 
The typical hopping length is then 



kg(fi)a 3 



P~ 1. 



f = a{T /T)^ 



(6.25) 



(6.26) 



(6.27) 



which is temperature dependent. This is why the conductance of the type discussed above 
is called the variable hopping conductance. The coefficient (3 can be found with the help of 
percolation theory and numerical calculations. For 3D case /3 f=a 21. 
For 2D case 

,(f /T)V* ™ P 



P = p e" 



P « 14. 



(6.28) 



kg 2 (p)a 2 

Here #2 is the two-dimensional density of states. Note that the derivation is based upon 
the assumption of a constant density of states at the Fermi level. We shall see that this 
assumption can be violated if electron-electron interaction is taken into account. 



6.5 Coulomb gap in the density of states 

At zero temperature, the distribution of electron is given by the minimum of electrostatic 
energy 



n 



1 d f nd ° n (l-n fc )(l-n fc Q _ d ° n ^ c (l-n fc 1 ™ c ^ 

k k'j^k kk 1 i ik i i'jti 



under condition that Ylk n k — const, or unconditional minimum of the functional 

H = H- pJ2 n >" 
k 



(6.29) 



(6.30) 



where the Fermi energy /1 should be then determined from the neutrality condition. We 
shall consider the energy counted from the chemical potential ii, 



6i = £i- p 



don (-1 acc 1 



rik 



j r ij 



-LI. 



(6.31) 



The goal is to find the distributions {n k } and {e^} that is a complicated many-body 
problem. We have discussed is previously for the limiting cases K <§C 1 and K ^> 1. Now 
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Density of states in a lightly doped 
and compensated semiconductor in the presence of a 
Coulomb gap at different degrees of compensation. 

\ the density of 
) *-=n.5, (b) 

K = t)A > {c)K = 0.9 



The energy £ is in units of e l Np 1 k ' 
states in N {e 2 N]? k V C (a 



Figure 6.2: 



we turn to the case K ~ 1 and the only energy parameter is eo = e 2 /M"D, where td is the 
average distance between the donors. The results of numerical simulations of the density 
of states for different K is shown in Fig. f>.'J\ . A gap at the Fermi level is clearly seen. The 
reason of this gap can be explained as follows. 

Consider a band (— e/2,+e/2) centered at the Fermi level. Let us estimate the energy 
cost of the electron transfer from the occupied state i to the empty one, j. Since the 
electron will be attracted to the donor, the energy cost is 



A? 



(6.32) 



It must be positive if the system was in the ground state. In other words, any two donors 
in this band having the energies in the opposite side of the Fermi level (e, ■ < 0) must 
be separated in space by a distance 



r ij > e 2 /ne . 



Consequently, the donor concentration, n(e) in a band of width e cannot exceed e 3 K 3 /e 6 . 
As a result, the density of states, g{l) = dn{e)/de must vanish when e — > 0, at least as fast 
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as e 



' 2 . Because of that, the usually accepted formula for the density of states is 



9\* 



') 



a 3 e 2 /t 3 /e 6 , d = 3 
a2\e\K 2 /e 4 , d = 2 



(6.33) 



Here ctj are numerical coefficients, d is the dimensionality (in 2D case the density of states 
is calculated per unit area). It is important that the density of states has a universal 
character. A sort of self-consistent treatment provides 0:2 = 2/71, o 3 = 3/n. 

The Coulomb gap modifies the variable range hopping conductance. In 3D case the 
Mott's law is modified as 




(6.34) 
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Chapter 7 

AC conductance due to localized 
states 



As we have seen, at low temperatures dc conductance vanishes in the insulating phase, the 
electronic states being localized. However, ac conductance remains finite. The reason is 
that an electron (or hole) can hop between the adjacent localized states giving rise to a 
finite energy dissipation. 

Deeply in the insulating phase the most probable is to find only pairs of localized states 
which are separated by the distances much less than the average. To analyze the physics 
of dissipation, let us consider such a pair. We used such a concept to discuss Anderson 
localization. 

Consider a close pair of localized states. As we have seen, the quantum states of an 
electron sharing the pair can be expressed through one-center wave functions as 

lj>± = Cl01 ± C 2 2 , 

which must satisfy matrix Schrodinger equation, 



A/2-E I 
I* -A/2 




0. 



Here the origin for energy is taken at (ei + 62)/% while A = e\ — e 2 . The secular equation 
is thus 

E 2 - (A/2) 2 - \I\ 2 = -> E± = ±y / (A/2) 2 + |J| 2 . 

Consequently, 

E + -E„ = W = V / A 2 + 4|/| 2 , (7.1) 

* = _?3eL. (7 . 2) 
c 2 A±W y ' 

Thus at A ^> I either c\ or c 2 is close to 1, and collectivization does not occur. 

To find the dissipation one has to calculate the contribution of a single pair, and then 
sum the contributions of different pairs. 
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Dissipation by an isolated pair 

One can discriminate between two absorption mechanisms. The first one is resonant ab- 
sorption due to direct absorption of photons accompanied by transition between the states 
tp + and t/>~. The second mechanism is due to phonon-assistant transitions - an external 
field modulates the occupation numbers of close pairs. The modulation lags in phase com- 
paring to the filed due to relaxation. As a result, a finite dissipation appears. Below we 
consider both mechanisms. 

Resonant contribution 

The resonant absorption is due to transition of an electron from the state with the energy 
E~ to the state E + . The energy absorbed by a pair per unit time due to an electric field 

£ = So cosut = -£ [exp(iut) + exp(— iut)] (7.3) 

can be written as (Fermi golden rule!): 

q= ^^| e £ (_| r | + )| 2 8 (to - W) (n- - n+) , (7.4) 



h 4 



Occupation numbers 



The occupation numbers n± are determined from the following considerations. Let us write 
down a two-site Hamiltonian 

e 2 

#1,2 = eifti + £2"-2 H n\n 2 + I(r)(ata 2 + at ax) . (7.5) 

nr 

Here rij are the occupation numbers, while q include Coulomb interaction with other 
neighbors. 



The Hamiltonian ( |7.5|) describes 4 states of the pair: 

1. The pair has no electrons. The energy E = 0. 

2. The pair has 1 electron. There are two states with energies 

1 2 2 

3. The pair has 2 electrons. There is one state with the energy 

e 2 

E 2 = ei + e 2 H • 

nr 
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Consequently, the probability to find a pair with 1 electron, the lower level to be occupied 
is 

i / £r-/A 

z = 1+exp (~^^J +exp r^^J +exp riff" 

The occupation number n + = n_ exp(— W/kT), and we obtain finally, 



q = YIz \ e£ o(-\^\+)\ 6 (ft" - w ) ex P — Jj^ 



x - exp (-^)j 



(7.6) 



(7.7) 



Relaxational contribution 

To analyze this contribution one has to consider the balance equation, say for n+(t), 

dn + n + — no(t) 
~df = t ' 
Here 

"°W = exp[^)V] + l - ^ = ^ [A + e£(i) ' r|2 + ^ ' (7 - 8) 
while r is the population relaxation time. Substituting 

n+(t) = n (t) + m(t) , ni(t) oc [e~ iujt + /i.e.] 

one easily obtains the relevant contribution to the absorbed energy as 

q = 7T dtW(t) ni (t). 
2tt Jo 

We get 

_ legp-if /A\ 2 u 2 r(W,r) 1 ( q 

q 2 \w) 1 + \uT{W,r)} 2 AkT cosh{W/2kT) ' { ' } 

The last factor is just — (dn /dW). 

To calculate the dissipation one has to specify the relaxation time, which depends in 
general on W and r. To do that let us specify the Hamiltonian to describe coupling 
between localized electrons and phonons. To construct the Hamiltonian, let us start with 
the unperturbed one, 

"•=sU> -aB- +/w - (7io) 

Here we introduce Pauli matrices, 

ai = ( 1 ' ° 2 = ( i j ' ^ = ( -1 ) ■ 
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Under influence of the phonon-induced strain the energy of each (j) component of the pair 
acquires the term proportional to the strain tensor, 



Thus, 



if duiiYj) du h {tj)\ 

«int = \ E [ A .* - ASWfo)] ff 3 , (7.11) 

where are the component of the deformational potential tensor for each component. 

Now we can make the transformation to a new basis which makes the Hamiltonian H,q 
diagonal. Technically, it is convenient to rewrite (|7. 1 1|) in the form 

_ , W f cos y sin v \ A 
^° = T v sin X -cos yJ' COSX =W 

Then one can represent the old basis through the new one as 

/ 0i \ ~ f ~ f cos x sin x \ * 

= T , T=\ . = 1 cos x - ^2 sm x • 

V 02 J J \ sm X cosx J 

Having in mind the algebra for Pauli matrices, 

a] = 1, a 2 cr 3 = ia u a 3 at = ia 2 , a x o 2 = ia 3 
we obtain the interaction Hamiltonian in a new basis, 



'hit ~ T~%nt f = \ E [ASW(rO - AL 2) « ife (r 2 )] ( A ff3 _ ^l^) . ( 7 .12) 



w w 



We are interested in the item proportional to o~\ which is responsible for phonon-assisted 
transitions between the levels. Using Fermi golden rule to calculate the relaxation rate, we 
get 

= f^V (7 13) 

Here we have extracted the coordinate-dependent factor (2I(r)/W) 2 = 2I exp(— 2r/a). 
The quantity r min has a transparent physical meaning of the minimal relaxation time for 
a pair with given interlevel spacing W. It is dependent on several characteristics of the 
electron-phonon interaction which we do not discuss now. 



Summation over the relevant pairs. 

Now we can proceed adding the contributions of different pairs that appears rather tricky. 
What we need is to find the so-called pair distribution function which is the probability to 
find a pair with bare energy spacing W and spatial distance r. 
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Figure 7.1: Energy scheme of absorption by a pair with the length r with allowance of 
Coulomb interaction. 



Resonant contribution 

For simplicity, let us assume low temperatures, T < fhu/k. Then n_ — n + « 1 at 

- fi < , - ji < E 2 - 2n . 

It means that one can introduce the variables A = ei — e 2 and £7f = (ei + e 2 — W)/2 
instead of ei and e 2 . One can show that V « 1 at 

e 2 

w < - // < o , 

and 

J dE^(n_ = n+) . . . -> + ^ . . . 

One can explain qualitatively this factor as follows from Fig. 0. When the center of gravity 
(ei + £2) /2 falls into the region of the width ~ e 2 / «r below the Fermi level the pair remains 
singly ionized and contains only 1 electron. 

The full scheme of calculations can be described as follows. We define 

a(uj) = 2 ^ 2 9t =jb-[dr[ de 2 g(e 2 ) f de x g{e x ) q(e x , e 2 , r) . (7.14) 

For the simplest (Anderson) model 

g(e)=g Q(A/2-\e\). 
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Another important simplification arises from exponential decay of the overlap integral, 
J(r) oc exp(— r/a). As a result, the coupling matrix element in the case of resonant 
contribution (as well as relaxation time in the case of the relaxational contribution) depend 
on r much stronger that other quantities. As a result, one can replace r apart from 
exponential functions by a proper characteristic length r c (which is different for the resonant 
and the relaxational contributions). 

Now we are ready to specify the results for low temperatures.. Substituting expression 



([7.1 5| ) for the contribution of a single pair we get (at T = 0) 



, , 87r 2 e 2 (7n f 00 A e 2 \ r 4 I 2 (r)dr 
a(u) = — — — / [hu+ 1 



3ft Jr u V Kr ) [(M 2 -4/ 2 (r)] 1/2 

2 / e 2 



l7[2 -e 2 glaurt( ktu + — ) (7.15) 



where 



r u = ahi(2Io/hu}) . (7.16) 
The case of finite temperatures is more complicated. At the region 

Tiuj kT e 2 1 'kTu 

the width of integration over E~ remains essentially the same. However the occupation 
numbers of the states Ef become very close, and 

(n_ - n + ) w tanh(^u;/2A;T) . 

This estimate is valid within all the region 

hto, kT <C e 2 / 'nr u . 

The case of very large temperatures can be considered in a similar way. At 

kT ^> e 2 / kt w , hu 

both states Ef fall into the layer < kT around the Fermi level, the difference (n_ — n + ) 
being 

1 

4cosh^[( J Ef -n)/2kT\ ' 



(n_ -n + ) ~ 2 



Thus. 
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Relaxational contribution 

Let us start with expression (|7.9p which assumes that there are only the pairs with one 
electron. Introducing the distribution function F(W, r) of the energy and spatial spacings 
(W and r) one obtains 

e 2 vr 2 [°° . F(W,r)dW r 

= 3*7^ A rdr Jo costf(W/2kT)l + ^- (7 - 17) 

We will see later that F(W, r) is a smooth function of r comparing to 

r(W,r) oc exp(2r/a) . 

Thus let us define 

r = t exp(2r/a) 

where tq is a smooth function of W and (in general) r. The properties of the ratio 

OJT utq exp(2r/a) 

1 + cj 2 r 2 1 + u; 2 Tq exp(4r/a) 

depend strongly on the product ut . At 

ujtq > 1 , r > a 

one has ^ 

G ~ exp(— 1r I a) . 

LUTq 

As a result, 

= 4^ f^dr f 1 . (7.18) 

V ; 3A;Tio io cosh 2 (W/2kT)T {W,r) v ! 

Thus, o appears u- independent, the temperature dependence being determined by the 
properties of the functions F(W, r) and tq(W, r). It is important that for the relevant pairs 
W ~ fcT, r = r T ~ (a/2) \n(2I /kT). 
At 

UOTq <^ 1 

the ratio G has a sharp maximum at 

r = r CUJ = (a/2) ln(l/wr ) 

and we can express G as 



In this way we obtain 



G « — 5(r - r c J . 

7T 



, , 4e 2 4 r F(iy,rjdiy 

ciom = — wor„. / s- ; oc . (7.19) 

v ; 3A;T ^io cosh 2 (W^/2A;T) V ; 
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Let us discuss a bit the distribution function F which is just the probability to find a single- 
electron pair with the energy distance W, spatial distance r. In the absence of Coulomb 
interaction, 

F(W, r) = J dei de 2 g(e 1 )g(e 2 ) 5 (w - y /(e 1 - e 2 ) 2 + 4/ 2 (r)) 
x {[/ (^r)[l " fo(Et)} + / (^ + )(l - MET)]} • 

Here /o is the Fermi function. The integration is very simple for the Anderson model 
(9 — 9o = const). In that case 

Coulomb interaction changes the situation, and one has to calculate the number of pairs 
taking into account electron repulsion at the same center. At low enough at kT <^ e 2 /nr cu} 
that leads to the result similar to the one discussed above, 

Thus at kT <^ e 2 /nr cu] the absorption is essentially temperature-independent at ut <^ 1. 



Part II 

Heavily Doped Semiconductors 
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Chapter 8 

Interband light absorption 



A typical experimental picture of frequency- dependent light absorption in heavily doped 
degenerate GaAs is shown in Fig. |8| The main features are: 

• rapid decrease of absorption (4 orders of magnitude) inside the forbidden gap (1.51 
eV at 77 K). 

• Shift towards short waves with increasing concentration (Moss-Burstein shift). 

The last feature is due to Pauli principle which requires that the transitions should take 
place only into empty states (above the Fermi level). For indirect transitions the shift is 
equal to fi while for direct ones it is /x(l + m e /mh), see Fig. [|. This picture is relevant to 
very low temperatures. At finite temperatures there is an absorption below the threshold 
because one can find a hole in the conduction band with the energy \x— e with the probability 
exp(— e/kT). Thus, 

a oc exp[— (// + EghuS) /kT] . 

Another important source which we are going to discuss are transitions from fluctuation 
levels above the top of valence gap, Fig. || We see that the absorption is proportional to 
the probability to find a level with the energy = E g + fi — huj which decays exponentially 
with €h- This probability must be multiplied by the transition probability which is not 
exponentially small in many important cases. 

Thus, the frequency dependence of absorption in degenerate doped semiconductors 
reproduces the DOS profile for minority carriers. 

The situation is somewhat different in the case of non-degenerate materials where the 
Fermi level is situated deep in the forbidden gap. The situation is realized at high temper- 
atures or degrees of compensation, as well as in indirect materials. 

Assume that there are impurities of both signs, N t = No + Na- The impurity potential 
is regarded as Coulombic at short distances r < r and screened at large distances. Let 
r be large enough to consider the potential as classical; we also neglect the correlation in 
defect positions. Mathematically, that can be expressed as 

h 2 h 2 p 2 
m e r^ m h r$ nr 
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Figure 8.1: The absorption coefficient as a function of the energy of light quanta at 77 K 
in n-GaAs. n ■ 1(T 17 , cm- 3 : 1 - 0.02, 2 - 2.2, 3 - 5.3, 4 - 12, 5 - 16.2, 6 - 31.5, 6 - 65. 
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Figure 8.2: Scheme of typical interband absorption processes. 
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Figure 8.3: Scheme of interband transitions in a degenerate semiconductor at T = 0. 
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Figure 8.4: Absorption of a quantum of deficit A in a non-degenerate semiconductor. 



One can show that in such the DOS tail can be represented as 

g{e) = 0(0) exp(-e 2 / 7 2 )- 

In such an important case the absorption behavior does not represent the behavior of 
DOS, see Fig. § Suppose we are studying the transition which results in the formation of 
an electron of energy e e and a hole of energy e^, 

A = e h - e e 

(here the energies are calculated from the non-perturbed band edges). We observe that the 
electron and hole are spatially separated. Thus the matrix element involves overlap integral 
which has almost nothing to do with DOS. What we meet here, is the Franz-Keldysh effect 
in a random field. 

To derive the result we use the conventional optimum fluctuation method. Here we 
demonstrate a simplified version of the method. 
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Consider a volume with linear dimension R. Charge fluctuations in the volume create 
a uniform electric field E, determined by the condition 

eER = A . 

The excess number Z of charge defects to create the field E = Ze 2 / ' kR 2 must be 

EkR 2 _ kRA 

Z — — — — . ( 

e e z 

The contribution of such a fluctuation to the absorption coefficient is proportional to 

1. the probability to find Z excess charges in the given volume, 

ex V (-Z 2 /N t R 2 )' , 

2. the probability for an electron to tunnel the distance R to meet the hole, 



5.1) 



exp y—R^Jm^A/hj 

(we assume m h 3> m 2 ). 
Substituting Z from Eq. (|3.1|) we obtain the probability to be proportional to 



exp 



Maximizing the exponent, we obtain 



- 2 A 2 R^^K" 



e 2 N t R 



n 



R = a 



Consequently, 



g(A) 
a(0) 



exp 



Ej (iV t a 3 )V2 • 
1 



where (3 is the number of the order 1. The derivation is valid at 

R<r . 

In the opposite case, the probability to find the proper fluctuation decreases, and the 
optimal cluster has the size r . In that situation one has to substitute r instead of R to 
obtain 

a(A) / A 2 \ 

exp 



a(0) 



T 
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